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The  Euclidean  or  Bunch-Davies  0(4, 1)  invariant  ‘vacuum’  state  of  quantum  fields  in  global 
de  Sitter  space  is  shown  to  be  unstable  to  small  perturbations,  even  for  a  massive  free  field 
with  no  self-interactions.  There  are  perturbations  of  this  state  with  arbitrarily  small  energy 
density  at  early  times  that  is  exponentially  blueshifted  in  the  contracting  phase  of  ‘eternal’ 
de  Sitter  space,  and  becomes  large  enough  to  disturb  the  classical  geometry  through  the 
semiclassical  Einstein  eqs.  at  later  times.  In  the  closely  analogous  case  of  a  constant,  uniform 
electric  field,  a  time  symmetric  state  equivalent  to  the  de  Sitter  invariant  one  is  constructed, 
which  is  also  not  a  stable  vacuum  state  under  perturbations.  The  role  of  a  quantum  anomaly 
in  the  growth  of  perturbations  and  symmetry  breaking  is  emphasized  in  both  cases.  In  de 
Sitter  space,  the  same  results  are  obtained  either  directly  from  the  renormalized  stress  tensor 
of  a  massive  scalar  field,  or  for  massless  conformal  fields  of  any  spin,  more  directly  from  the 
effective  action  and  stress  tensor  associated  with  the  conformal  trace  anomaly.  The  anomaly 
stress  tensor  shows  that  states  invariant  under  the  0(4)  subgroup  of  the  de  Sitter  group 
are  also  unstable  to  perturbations  of  lower  spatial  symmetry,  implying  that  both  the  0(4, 1) 
isometry  group  and  its  0(4)  subgroup  are  broken  by  quantum  fluctuations.  Consequences 
of  this  result  for  cosmology  and  the  problem  of  vacuum  energy  are  discussed. 
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I.  ‘VACUUM’  STATES  IN  DE  SITTER  SPACE 


The  existence  of  a  ground  state  as  the  state  of  lowest  energy  is  fundamental  to  all  quantum  mechanical 
systems.  For  quantum  field  theory  (QFT)  in  flat  Minkowski  spacetime  the  vacuum  state  is  defined  as 
the  eigenstate  of  the  Hamiltonian  operator  of  the  system  with  the  lowest  eigenvalue.  The  existence  of  a 
Hamiltonian  generator  of  time  translational  symmetry,  with  a  non-negative  eigenvalue  spectrum,  bounded 
from  below  is  crucial  to  the  existence  and  determination  of  the  vacuum  ground  state. 

This  definition  of  the  vacuum  in  flat  spacetime  makes  use  of  an  essential  property  of  the  Poincare 
group,  namely  that  positive  and  negative  (particle  and  antiparticle)  halves  of  the  Hamiltonian  spectrum 
do  not  mix,  remaining  distinct  under  any  of  the  continuous  generators  of  the  group.  Hence  the  vacuum 
state  in  flat  space  QFT  is  the  same  for  all  inertial  frames  related  to  each  other  by  translations,  rotations 
and  Lorentz  boosts,  and  the  vacuum  enjoys  complete  invariance  under  Poincare  symmetry. 

As  is  well  known,  none  of  these  properties  hold  in  a  general  curved  spacetime,  in  time  dependent 
background  fields,  nor  even  in  flat  spacetime  under  general  coordinate  transformations  which  are  not 
Poincare  symmetries.  In  these  circumstances  the  definitions  of  ‘vacuum’  and  ‘particles’  become  much 
more  subtle.  Related  to  this,  whereas  the  infinite  zero  point  energy  associated  with  the  QFT  vacuum 
may  be  disregarded  as  unobservable  in  flat  space  QFT,  the  energy  of  the  quantum  vacuum  in  curved 
spacetime  cannot  be  neglected  when  the  coupling  to  gravity  is  taken  into  account. 

These  issues  come  to  the  fore  in  the  important  special  case  of  de  Sitter  space,  the  classical  spacetime 
with  a  positive  cosmological  constant  A  >  0,  which  itself  may  be  regarded  as  a  vacuum  energy  density 
uniformly  curving  space.  The  geodesically  complete  full  de  Sitter  manifold  may  be  represented  as  a  single 
sheeted  hyperboloid  of  revolution  embedded  in  five  dimensional  flat  Minkowski  spacetime,  c.f.  Fig. 

[T].  It  has  the  isometry  group  0(4, 1)  with  10  continuous  symmetry  generators,  the  same  number  as  the 
Poincare  group  of  Minkowski  space,  and  the  maximal  number  possible  for  any  solution  of  the  vacuum 
Einstein  field  equations  in  4  spacetime  dimensions.  This  maximal  symmetry  is  evident  from  the  constant 
and  uniform  Riemann  and  Ricci  curvature  tensors  and  scalar  of  de  Sitter  space,  which  are  respectively 


(1.1a) 

(1.1b) 

R  =  =  4A . 

(1.1c) 

A  natural  attempt  to  generalize  the  QFT  vacuum  of  flat  space  to  de  Sitter  space  makes  use  of  this 
geometrical  symmetry  of  de  Sitter  space  to  define  the  de  Sitter  invariant  ‘vacuum’  |t;)  as  the  state 
possessing  the  same  maximal  0(4, 1)  symmetry  in  the  Hilbert  space  of  states.  Introduced  by  Chernikov 
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FIG.  1.  The  de  Sitter  manifold  represented  as  a  single  sheeted  hyperboloid  of  revolution  about  the  T  =  axis, 
embedded  in  five  dimensional  flat  spacetime  X°‘),  a  =  1, ...  4,  in  which  the  X^ ,  X^  coordinates  are  suppressed. 
The  hypersurfaces  at  constant  T  =  X^  =  sinhu  are  three-spheres,  S^.  The  at  T  =  ±oo  are  denoted  by  I±. 

and  Tagirov  (CT)  [2],  this  state  is  commonly  known  also  as  the  Bunch-Davies  (BD)  state  OH],  or  the 
Euclidean  ‘vacuum,’  because  its  Green’s  functions  are  those  obtained  by  analytic  continuation  from  the 
Euclidean  at  least  for  massive  fields  where  no  obvious  infrared  issues  arise  [5]. 

It  is  important  to  recognize  that  unlike  in  flat  space,  the  construction  of  the  CTBD  state  is  not  based 
on  diagonalization  of  any  Hamiltonian  nor  any  minimization  of  energy.  In  fact  no  suitable  Hamiltonian 
operator  with  a  spectrum  bounded  from  below  exists  at  all  in  de  Sitter  space,  even  for  free  QFT.  In  the 
globally  complete  coordinates  of  the  de  Sitter  hyperboloid 

ds'^  =  [—du^  +  cosh^udS^)  (1.2) 

with 

=  dN  ■  dN  =  dx^  +  sin^x  dfi  ■  dn  =  dx^  +  sin^x  +  sin^d  d0^)  (1.3) 

the  standard  round  metric  on  S^,  the  de  Sitter  metric  is  dependent  on  the  time  u.  Thus  translation  in 
u  is  not  a  symmetry  of  de  Sitter  space  and  the  generator  of  u  time  translations  is  not  conserved.  As  a 
consequence,  the  ‘vacuum’  defined  by  Hamiltonian  diagonalization  at  one  instant  of  u  time  will  contain 
‘particles’  at  any  other  u  time.  This  is  equally  true  in  the  flat  spatial  slicing  of  de  Sitter  space 

ds^  =  — dr^  -b  dx  ■  dx  (1.4) 

used  most  frequently  in  cosmology,  which  is  similarly  dependent  on  the  time  r. 

The  non-existence  of  a  conserved  Hamiltonian  generator  bounded  from  below  in  de  Sitter  space  is 
a  consequence  of  the  de  Sitter  symmetry  group  0(4, 1)  itself.  Unlike  the  Poincare  group,  no  invariant 
separation  into  positive  and  negative,  particle  and  antiparticle  subspaces  exists  in  de  Sitter  space,  and 
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any  de  Sitter  symmetry  generator  chosen  for  the  role  of  the  Hamiltonian  has  a  spectrum  of  both  positive 
and  negative  eigenvalues  which  are  mixed  by  the  action  of  other  generators  of  the  group  [6] .  One  of  the 
4  non-compact  Lorentz  boost  generators  of  the  0(4, 1)  symmetry  group  may  be  selected  (arbitrarily)  as 
the  Hamiltonian  of  the  system,  generating  time  translations  t  — )•  f  -|-  At  in  the  static  coordinates  of  de 
Sitter  space,  where  the  line  element  takes  the  form 

ds^  =  —(1  —  dt^  H - „  -|-  r^{d6‘^  +  siT?9  d4>^) .  (1.5) 

1  — 


In  these  coordinates  the  geometry  is  independent  of  the  time  t.  However  the  event  horizon  at  r  =  H~^ 
relative  to  the  origin  r  =  0  is  now  manifest,  and  the  static  coordinates  cover  only  one  quarter  of  the 
full  de  Sitter  manifold.  The  Killing  symmetry  d/dt  is  not  globally  timelike,  and  changes  its  orientation 
from  one  quadrant  to  another,  as  may  be  seen  from  the  Carter-Penrose  conformal  diagram  of  de  Sitter 
space:  Fig.  A  direct  consequence  of  this  is  that  the  corresponding  Hamiltonian  symmetry  generator 
across  any  complete  Cauchy  surface  is  not  positive  definite,  but  rather  unbounded  from  below,  as  Lorentz 
boosts  are.  Hence  its  eigenstates  or  expectation  values  cannot  be  used  to  select  a  global  minimum  energy 
vacuum  state.  The  choice  of  d/dt  is  also  arbitrary  and  the  separation  into  positive  and  negative  energies 
with  respect  to  d/dt  is  non-invariant  under  de  Sitter  group  transformations.  The  particle  concept  is 
likewise  affected,  as  the  CTBD  de  Sitter  invariant  ‘vacuum’  state  |u)  is  actually  a  state  with  a  thermal 


distribution  of  ‘particles’  with  respect  to  the  Killing  Hamiltonian  generator  dt  of  (1.5)  with  the  Hawking 
de  Sitter  temperature  [7] 

hH 


T  = 

"  2TTk, 


(1.6) 


and  in  that  sense  is  not  a  vacuum  state  at  all. 

The  horizon  and  causal  structure  of  de  Sitter  space  raises  questions  of  how  a  vacuum  state  can  be 
prepared  operationally  even  in  principle.  Inspection  of  the  conformal  diagram  in  Fig. shows  that  points 
on  a  Cauchy  surface  at  uq  <  0  with  widely  different  N  (for  example  at  y  =  0  and  its  antipodal  point 
X  =  tt)  could  never  have  been  linked  by  any  causal  signal  in  the  past.  As  the  initial  time  uq  is  taken  earlier 
and  earlier,  this  causal  disconnection  affects  more  and  more  of  the  initial  u  =  uq  Cauchy  surface.  Since 
past  inhnity  /_  is  spacelike,  as  uq  — )>  — oo  in  this  limit  no  two  different  points  on  could  have  been  in 
any  causal  contact  whatsoever.  Thus  any  global  initial  data  on  §^,  including  that  necessary  to  construct 
the  CTBD  state  |u)  cannot  have  been  provided  at  any  initial  time  uq  <  0  by  any  causal  process  within  de 
Sitter  space  itself.  Instead  initial  data  has  simply  to  be  posited  over  the  full  spacelike  S^,  at  points  outside 
of  the  causal  horizon  of  any  local  agent  who  might  have  prepared  it  at  early  times.  This  is  equivalent  to 
the  existence  of  a  particle  horizon  [T],  and  is  completely  unlike  that  of  flat  Minkowski  spacetime,  where 
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FIG.  2.  The  Carter-Penrose  conformal  diagram  for  de  Sitter  space,  in  which  light  rays  emanating  from  any  point 
are  at  45°,  and  the  angular  coordinates  9,(j)  are  suppressed.  The  quarter  of  the  diagram  labeled  as  the  static 


region  are  covered  by  the  static  coordinates  of  (1.5).  The  orbits  of  the  static  time  Killing  field  d/dt,  r  =  0 


and  curves  of  constant  r  >  0  are  shown.  In  contrast,  the  surfaces  of  constant  u  time  coordinate  of  (1.2)  are 
horizontal  straight  lines  across  the  diagram  with  x  G  [0,7r],  labelled  on  the  right  by  the  conformal  time  coordinate 
r]  =  sin“ ^  (tanh  u)  S  (— f ,  f )  as  u  ranges  from  — oo  to  +oo  at  past  (/_)  and  future  (/+)  infinity  respectively. 


Cauchy  data  on  a  fixed  time  slice  t  =  const,  may  be  prepared  in  principle  by  transmitting  signals  causally 
from  a  single  point  early  enough  in  the  past.  The  existence  of  horizons  and  the  absence  of  a  global  time 
coordinate  connected  with  any  symmetry  reveal  the  essential  difficulties  with  defining  a  global  vacuum 
state  for  QFT  in  eternal  de  Sitter  space. 

They  also  imply  that  the  mathematical  requirement  of  global  de  Sitter  invariance  cannot  be  realized 
by  any  local  physics  within  de  Sitter  space  itself,  requiring  instead  an  acausal  fine  tuning  of  initial  data 
at  spacelike  past  inhnity  with  a  view  to  the  entire  future  manifold,  which  is  presumed  to  be  known 
in  advance  in  order  to  specify  a  globally  0(4, 1)  invariant  state.  Although  maximal  0(4, 1)  symmetry 
may  seem  natural  mathematically,  or  by  analytic  continuation  from  the  Euclidean  where  no  causal 
relations  apply,  it  is  quite  unnatural  with  respect  to  the  physical  principles  of  locality  and  causality  in 
real  time,  as  well  as  lacking  any  de  Sitter  invariant  Hamiltonian  minimization  principle. 

For  these  reasons  it  is  important  to  study  the  sensitivity  of  physical  quantities  in  de  Sitter  space 
to  fluctuations  and/or  perturbations  of  states  away  from  precisely  the  ‘right’  one  for  global  0(4,1) 
invariance,  rather  than  simply  assuming  this  symmetry.  The  calculation  of  the  imaginary  part  of  the 
effective  action  of  a  simple  scalar  QFT  in  de  Sitter  space  due  to  particle  creation  urn  already  shows 
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that  de  Sitter  space  is  unstable  to  spontaneous  creation  of  particle  pairs  from  the  vacuum,  just  as  is  an 
‘eternal’  uniform  electric  field  'Ei  =  Ez,  permeating  all  of  space  m-  This  electric  field  analogy  and  the 
close  relation  between  fluctuations  and  dissipation  in  any  causal  theory  suggests  that  the  ‘shorting  of  the 
vacuum’  should  result  in  the  classical  energy  of  de  Sitter  space  converting  itself  into  standard  matter  and 
radiation,  thus  providing  a  route  to  a  dynamical  solution  to  the  cosmological  ‘constant’  problem  [mill]. 

As  in  electrodynamics,  interactions  in  de  Sitter  space  are  certainly  relevant  to  understanding  of  the 
detailed  evolution  and  final  state,  particularly  since  spontaneous  pair  creation  should  be  accompanied 
by  induced  emission  processes  which  can  create  an  avalanche  of  particles  that  will  inevitably  interact 
and  thermalize,  leading  to  the  final  dissipation  of  vacnum  energy  into  matter  and  radiation.  A  fuller 
understanding  of  these  non-equilibrium  processes  may  well  lead  to  a  satisfactory  resolution  of  the  cosmo¬ 
logical  ‘constant’  problem,  and  be  relevant  to  observational  cosmology  through  the  residual  dark  energy 
in  the  present  epoch  mm-  However,  this  dynamics  has  not  been  fnlly  solved  in  in  four  dimensions 
even  in  flat  space  electrodynamics.  Moreover  a  number  of  questions  persist  about  QFT  in  de  Sitter 
space,  even  in  the  non-interacting  case,  and  these  should  be  settled  definitively  first,  because  they  depend 
through  the  energy-momentum-stress  tensor  Tab  only  npon  the  universal  coupling  to  the  gravitational 
field,  independently  of  any  matter  self-interactions. 

In  this  paper  we  study  the  behavior  of  the  renormalized  energy-momentum  tensor  {Tab)  of  QFT  under 
perturbations  of  the  |u)  state  to  nearby  states  of  lower  symmetry.  In  the  expanding  part  of  de  Sitter  space 


n  >  0  of  ( 1.2 )  or  in  the  Poincare  coordinates  of  flat  spatial  sections  ( 1.4 ),  it  has  been  shown  that  in  a  fixed 
de  Sitter  background,  the  expectation  value  {Tab)  for  a  scalar  field  with  effective  mass  +  ^R  >  0 

approaches  an  0(4, 1)  de  Sitter  invariant  value  at  late  times,  for  all  spatially  homogeneous  UV  allowed 
perturbations  m-  Physically  this  result  may  seem  intuitively  obvious,  since  all  deviations  from  the 
expectation  value  in  de  Sitter  invariant  state  |u)  are  redshifted  in  the  de  Sitter  expansion  and  vanish  in 
the  n  — )•  oo  limit.  Since  global  or  ‘eternal’  de  Sitter  space  is  time  reversal  invariant,  the  attractor  behavior 
in  the  expanding  phase  implies  just  the  opposite  behavior  under  time  reversal  in  the  contracting  phase. 
That  is,  very  small  changes  in  the  initial  state  in  the  very  distant  past  uq  — )>  — oo  of  eternal  de  Sitter 
space  with  initially  very  small  {Tab)  must  necessarily  produce  larger  and  larger  effects  in  {Tab)  as  the 
contraction  proceeds  towards  u  =  0.  This  is  just  the  case  where  the  aforementioned  issues  with  causality 
at  spacelike  /_  arise,  and  this  sensitivity  to  initial  conditions  at  /_  is  the  source  of  the  instability. 

By  studying  the  general  behavior  of  the  renormalized  {Tab)  in  states  with  lower  symmetry,  we  show 
in  this  paper  that  the  CTBD  de  Sitter  invariant  state  |u)  is  unstable,  in  the  sense  that  there  is  a  large 
class  of  initial  state  perturbations  which  have  exponentially  small  energy  density  in  the  infinite  past 
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uq  — >■  — oo  but  which  grow  large  enough  through  exponential  blueshifting  proportional  to  a~^,  where 
a  =  coshu  is  the  scale  factor  in  (1.2),  to  exceed  the  classical  background  energy  A/SttG  and  hence 
significantly  disturb  the  de  Sitter  geometry  at  tt  =  0.  In  fact,  there  are  such  states  with  (Tab)  larger  than 
any  finite  value  at  u  =  0. 

This  extreme  sensitivity  to  initial  conditions  as  uq  — >■  — oo  implies  that  0(4, 1)  de  Sitter  invariance  is 
broken,  and  the  spacetime  will  generally  depart  from  de  Sitter  space  when  the  backreaction  of  {Tab)  of 
any  matter  or  radiation  on  the  geometry  is  taken  into  account,  through  the  semiclassical  Einstein  eqs.. 


i2“,-|n  +  An  =  8vrO(n), 


(1.7) 


and  quite  apart  from  any  matter  self-interactions  or  higher  loop  effects.  Although  in  a  fixed  de  Sitter 
background  the  energy  density  of  spatially  homogeneous  perturbations  will  begin  to  decrease  again  for 


u  >  0,  perturbations  of  the  CTBD  state  and  their  backreaction  through  (1.7)  will  have  already  drastically 
altered  the  geometry  in  the  contracting  phase  and  broken  the  de  Sitter  symmetry  by  n  =  0,  rendering 
further  evolution  ignoring  backreaction  moot.  This  large  backreaction  of  the  energy- momentum  tensor 
for  perturbations  of  the  CTBD  state  is  independent  of  any  definition  of  particles. 

Although  for  definiteness  we  study  this  growth  of  (Tab)  explicitly  in  a  scalar  field  theory,  the  result 
is  clearly  much  more  general.  A  very  useful  tool  for  characterizing  the  behavior  of  the  stress  tensor  in 
any  coordinates  is  the  one-loop  effective  action  of  the  trace  anomaly  and  the  stress  tensor  derived  from 


it  |14H16j.  The  non-local  form  of  this  effective  action,  c.f.  (5.1)  already  indicates  infrared  de  Sitter 
breaking  effects,  and  sensitivity  to  initial  and/or  boundary  conditions  for  conformal  QFT’s  of  any  spin. 
The  corresponding  stress  tensor  may  be  found  in  closed  form  in  de  Sitter  space  in  any  coordinates  by 


solving  a  classical,  linear  eq.  (5.6)  for  a  scalar  condensate  effective  field,  whose  solutions  necessarily 
break  de  Sitter  invariance,  and  allow  wide  classes  of  initial  state  perturbations  for  fields  of  any  spin  to 
be  surveyed  at  once.  Because  de  Sitter  space  is  conformally  flat,  this  anomaly  stress  tensor  is  a  complete 
description  of  the  full  QFT  stress  tensor  for  conformal  fields  linearized  around  the  CTBD  state  |u)  at  all 
length  scales  much  larger  than  the  Planck  length  Lpi,  where  semiclassical  methods  should  apply  HZj. 

The  blueshifting  of  the  energy  density  of  even  massive  fields  to  eventually  ultrarelativistic  behavior 
shows  that  the  conformal  anomaly  stress  tensor  is  relevant  for  long  time  evolution  even  if  the  underlying 
QFT  is  not  conformally  invariant.  When  the  scalar  perturbations  are  spatially  inhomogeneous  new  effects 
may  also  become  apparent.  In  Ref.  m  we  studied  spatially  inhomgeneous  scalar  perturbations  in  linear 
response  of  conformal  QFT’s  around  de  Sitter  space  and  found  a  class  of  gauge  invariant  perturbations, 
which  do  not  redshift  away  but  instead  give  diverging  energy-momentum  components  at  r  =  H~^  in  static 
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coordinates  (1.5).  These  may  be  interpreted  as  fluctuations  in  the  Hawking  de  Sitter  temperature  (1.6) 


at  the  de  Sitter  horizon  with  respect  to  some  arbitrary  but  fixed  choice  of  origin,  and  clearly  respect  only 
rotational  0(3)  invariance  around  r  =  0  and  static  time  t  translational  invariance.  This  result  suggests 
that  fluctuations  on  the  horizon  scale  H~^  may  produce  significant  backreaction  in  de  Sitter  space,  and 
that  the  0(4, 1)  symmetry  is  unstable  to  such  spatially  inhomogeneous  scalar  fluctuations  in  the  Hawking 
de  Sitter  temperature  [18] .  Tensor  perturbations  have  been  studied  recently  in  m- 


Using  the  anomaly  form  of  Tab  we  shall  show  that  there  is  even  greater  sensitivity  to  spatially  inho¬ 
mogeneous  non-0(4)  invariant  initial  data  in  the  distant  past  uq  — >■  — oo  of  global  de  Sitter  space,  so  that 
0(4)  invariance  is  broken  as  well  as  the  full  0(4, 1)  de  Sitter  invariance.  This  strongly  suggests  that  spa¬ 
tial  inhomogeneities  are  more  important  in  QFT  in  de  Sitter  space  than  previously  suspected,  supporting 
the  results  of  m  Such  spatially  inhomogeneous  perturbations  clearly  are  relevant  even  in  the  expanding 


Poincare  patch  (1.4).  The  interesting  questions  of  the  behavior  of  the  stress  tensor  in  states  of  lower 


symmetry,  such  as  the  0(3)  symmetry  evident  in  static  coordinates  (1.5),  and  consequences  for  spatially 
inhomogeneous  cosmologies  will  be  taken  up  in  future  publications.  An  accompanying  and  closely  related 
paper  gives  a  fuller  treatment  of  the  instability  of  global  de  Sitter  space  to  particle  creation  [9| . 


The  paper  is  organized  as  follows.  In  the  next  section  we  construct  the  time  symmetric  invariant 
state  analogous  to  the  CTBD  state  in  de  Sitter  space,  in  the  case  of  a  uniform,  constant  electric  field 
background  E  =  i?  z,  and  show  that  it  also  is  unstable  to  perturbations  for  which  the  mean  current  {jz) 
grows  with  time.  This  growth  of  the  current  and  breaking  of  background  symmetries  can  be  understood 
by  consideration  of  a  quantum  anomaly,  in  this  case  the  chiral  anomaly  of  massless  fields  in  two  spacetime 
dimensions.  The  reader  interested  primarily  in  de  Sitter  space  proper  may  skip  this  section  upon  first 


reading  and  proceed  directly  to  Sec.  HI  where  we  begin  discussion  of  the  CTBD  state  and  general  states 


of  0(4)  symmetry  in  de  Sitter  space.  In  Sec.  IV  we  construct  the  renormalized  expectation  value  of 
the  stress  tensor  of  a  massive  scalar  field  with  conformal  coupling  ^  =  g  in  general  0(4)  invariant  states 
in  the  global  hyperboloid  coordinates  0  of  de  Sitter  space,  and  explicitly  exhibit  the  class  of  states 
with  large  backreaction  ai  u  =  0.  In  Sec.  |V]we  consider  the  effective  action  and  stress  tensor  associated 
with  the  trace  anomaly  of  conformal  fields  in  de  Sitter  space  and  show  how  the  strong  infrared  effects, 
sensitivity  to  initial  conditions,  and  breaking  of  de  Sitter  symmetry  is  inherent  in  the  conformal  anomaly 
for  QFT’s  of  any  spin.  In  Sec.  |VI|  we  extend  the  analysis  of  the  anomaly  stress  tensor  to  states  of  lower 
than  0(4)  symmetry,  showing  that  these  spatially  inhomogeneous  perturbations  grow  even  more  rapidly 


to  larger  values  at  u  =  0  than  0(4)  symmetric  states.  Sec.  VH  contains  our  conclusions  and  a  discussion 


of  their  possible  consequences  for  cosmology  and  the  problem  of  cosmological  vacuum  energy. 


II.  CONSTANT  UNIFORM  ELECTRIC  FIELD:  INVARIANT  STATE  AND  INSTABILITY 


A.  Time  Symmetric  Invariant  State 


The  example  of  a  charged  quantum  field  in  the  background  of  a  constant  uniform  electric  field  has 
many  similarities  with  the  de  Sitter  case.  Although  this  problem  has  been  considered  by  many  authors 
[TOl  [20l-[26]  ■  the  existence  of  a  time  symmetric  state  analogous  to  the  CTBD  state  in  de  Sitter  space  does 
not  appear  to  have  received  previous  attention,  and  is  particularly  relevant  to  our  study  of  vacuum  states 
in  de  Sitter  space,  so  we  consider  this  case  first  in  some  detail. 


Analogous  to  choosing  global  time  dependent  coordinates  (1.2)  in  de  Sitter  space,  one  may  choose  the 
time  dependent  gauge 


Az  —  —Et , 


At  —  Afr,  —  An,  —  0 


(2.1) 


in  which  to  describe  a  fixed  constant  and  uniform  electric  field  E  =  i?  z  in  the  z  direction.  Treating  the 
electric  field  as  a  classical  background  field  analogous  to  the  classical  gravitational  field  of  de  Sitter  space, 
the  wave  equation  of  a  non-self-interacting  complex  scalar  field  $  is 


[—{dfj,  —  ieA^){d^  —  ieA^)  +  m^]  $  =  0 


(2.2) 


in  the  classical  electromagnetic  potential  (2.1). 


The  solutions  of  (2.2)  may  be  decomposed  into  Fourier  modes  $  ~  ■^/k(t)  with 


^  /k(i)  =  0 


where  the  frequency  function  is  defined  by 


=  [{kz  +  eEt)"^  +  k\+  m^]  ^  =  V2eE  y  —  +  X. 


(2.3) 


(2.4) 


We  have  defined  the  dimensionless  variables 


u  = 


{kz  +  eEt) 


^  „  >  0 


(2.5) 


eE  '  ^  2eE 

and  chosen  the  sign  of  eE  to  be  positive  without  loss  of  generality.  With  /k(^)  — >•  fx{u),  the  dimensionless 


mode  eq.  (2.3)  becomes 


d?  V? 


f\{u)  =  0 


(2.6) 


the  solutions  of  which  may  be  expressed  in  terms  of  confluent  hypergeometric  functions  iTi(a;c;2;)  or 
parabolic  cylinder  functions  Dy  m- 
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Since  (2.6)  is  real  and  symmetric  under  u  — >■  —u,  it  is  clear  that  its  real  solutions  can  be  classified 


into  those  which  are  either  even  and  odd  under  this  discrete  reflection  symmetry.  Let  us  define  two 


fundamental  real  solutions  of  (2.6)  f^\u),i  =  0, 1  by  the  conditions 


[friu)r  =  friu),  z  =  o,i 

(real) 

(2.7a) 

II 

o 

T 

(even) 

(2.7b) 

II 

1 

T 

(odd) , 

(2.7c) 

which  are  even  or  odd  respectively,  and  which  satisfy  the  initial  data 

/f(0)  =  l,  /®'(0)  =  0, 

/f)(0)  =  0,  /l'^'(0)  =  l 

at  n  =  0,  where  the  primes  denote  differentiation  with  respect  to  u.  These  fundamental  real  solutions  of 


(2.8a) 

(2.8b) 


(2.6)  are  most  concisely  expressed  in  terms  of  the  confluent  hypergeometric  (Kummer)  function 


,  ,  .  ^  (a)n  z"'  ,  .  r(a  +  n) 

^>(a,  c;  z)  =  iFi(a;  c;z)  =  J2  TTT  TT  >  = 


n=o 


which  has  the  integral  representation  m 

r(c)  '■i 


^(a,c;z)  = 


r(a)r(c-  a)  Jo 


dxe^^  (1  -x) 


c—a—1 


Via) 


Re  c  >  Re  o  >  0 


(2.9) 


(2.10) 


in  the  form 


„iui/  X  ^  I  ^  iX  1  iu^\  ijd  / 1  zA  1  iu^\ 

=  e- 4  =e4  $  , 


do), 

fx\u)  =ue-"f  $  (  ^ 


(2.11a) 
(2.11b) 

are  clearly  even  or  odd  respectively,  real  by  the  Kummer  transformation  which  yields  the  second  forms 


iX 

3 

iu'^X 

iX 

3 

— 

=  u  e  4  cp  — 

2  ’ 

2  ■ 

2' 

2  y 

V4 

2’  2 

in  (2.11),  and  satisfy  the  initial  data  (2.8) 


(■i) 

It  is  also  possible  to  express  these  fundamental  real  solutions  as  linear  combinations  of  parabolic 
cylinder  functions  in  the  forms  m 



A-!r(f  +  t) 


vr 


1  iA(e  u)  +  D_i_^J-e  4  n) 


f[^\u)  =  2'2  ie  4 


„r(i  +  f) 


TT 


—  \  /  — 
-D-i-ixie  4  u)  +  D_i_.J-e  4  u) 


(2.12a) 

(2.12b) 


which  representations  are  useful  for  identifying  their  relationship  to  the  in  and  out  positive  frequency 
scattering  solutions  defined  as  u  ^  ^oo  respectively  in 
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From  the  fundamental  real  solutions  ( 2.8 )-( 2.11)  one  can  construct  the  complex  mode  functions 


vx{u)  =  {8eEX)  4  f^°\u)  -  iX^  f^\u) 


1  r\  ry  1  ivP' 

=  2“2  +  m  )“4  e~~ 


^  ,  1  zA  1  iu^ 
^  (  4  +  Y’  2’  ^ 


1  ,  /3  zA  3  iu^ 


which  are  normalized  according  to  the  Wronskian  condition 



and  which  satisfy  the  time  reversal  conjugation  property 


(2.13) 


(2.14) 



(2.15) 


These  v\  mode  functions  satisfy  the  initial  data 

1 


^^a(0)  =  2  2{kj_+m‘^)  4  = 


•\/2cuk 


u=0 


dv\ 

dt 


= - ^  {k\  + = -iu}k.vx{u)  (2.16) 

u=o  x/2 


which  coincides  with  the  definition  of  the  lowest  order  adiabatic  frequency  mode  functions  at  the  symmet¬ 


ric  point  zz  =  0.  The  solution  of  (2.6)  satisfying  conditions  (2.14)-(2.16)  is  unique.  Because  of  relations 


(2.12)  and  the  simple  asymptotic  forms  of  the  functions,  the  symmetric  mode  function  vx  is  a  coherent 


superposition  of  positive  and  negative  frequency  (particle  and  anti-particle)  solutions  as  zz  — Too,  just 
as  the  CTBD  mode  function  is  in  de  Sitter  space  |9]. 


The  existence  of  such  a  time  reversal  invariant  solution  to  ( 2.6 )  is  related  to  the  existence  of  a  maximally 


symmetric  state  constructed  along  the  lines  of  the  maximally  0(4, 1)  invariant  CTBD  state  in  the  de  Sitter 
background.  If  the  charged  quantum  field  $  is  expanded  in  terms  of  these  symmetric  basis  functions  in 
a  finite  volume  V 


1 



—  ikx 


(2.17) 


with  u  and  A  defined  by  (2.5),  then  a  symmetric  state  |z;)  in  a  background  constant  uniform  electric  field 
may  be  defined  by 


aj(|z;)  =  6K|z;)  =  0  .  (2-18) 

The  symmetry  in  this  case  is  isomorphic  to  the  full  Poincare  symmetry  group  of  zero  electric  field  in 
flat  Minkowski  space.  This  is  due  to  the  remarkable  fact  that  a  canonical  transformation  exists  that 
transforms  the  algebra  of  position  and  momentum  operators,  and  Pi^  =  —idjdx’^,  in  a  constant. 
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uniform  E  field  background  to  new  position  and  momentum  operators,  and  Pu  =  —idjdX'^,  such 


that  the  Klein-Gordon  operator  (2.2) 


-{dfj,  -  ieA^){d^  -  ieA^)  +  w?  =  -pf  +  {p^  +  eEt)"^  +  pI  +  pI  + 

=  -P|  +  P|  +  Pi  +  +  m2  =  0 


(2.19) 


(with  Px  =  Px,  Py  =  Py  and  Pz  =  Pz)  becomes  that  of  flat  space  with  zero  field  [28].  The  existence  of  this 
transformation  and  symmetry  may  be  less  surprising  when  it  is  recognized  that  there  are  two  quantities 


Pt  =  {Pt  -  ‘^eEt  —  e^E'^t^)  2  =  {pi  +  pi  +  p^  +  m^)  2 

TPz  +  ZPt  =  [eEz  +  pt  —  Pt) 
eE 


(2.20a) 

(2.20b) 


that  are  conserved  by  virtue  of  the  eq.  of  motion  (2.19),  and  (together  with  Pz  =  Pz  which  generates 


space  Z  translations)  they  generate  T  time  translations  and  Lorentz  boosts  in  the  Z  direction  in  the 
transformed  {T,Z)  coordinates.  This  dynamical  maximal  Poincare  symmetry  in  the  constant,  uniform 
E  field  is  analogous  to  the  maximal  0(4, 1)  point  symmetry  group  of  de  Sitter  space.  In  each  case  the 
existence  of  a  maximally  symmetric  state  \v)  which  enjoys  the  full  symmetries  of  the  background  follows. 

The  expectation  value  of  the  electric  current  operator  is  given  in  the  symmetric  state  |u)  by 


{'u\jz\v)n  =  2e 


d^k 


{kz  +  eEt) 


kA(w)P  - 


1 


2a;k(t)_ 


(2.21) 


j  (2vr)^ 

where  the  second  term  is  the  lowest  order  adiabatic  vacuum  subtraction  sufficient  for  the  constant  E 
field  background  |26].  Actually  by  changing  integration  variables  from  kz  to  u  and  using  the  fact  that 


both  |uA(M)r  and  Wk(t)  are  even  functions  of  u,  it  is  clear  that  both  terms  in  the  integrand  of  (2.21)  are 


odd  under  u  — )•  —u  and  thus  give  vanishing  contributions  if  integrated  symmetrically  in  u.  Hence  as  a 
consequence  of  time  reversal  invariance  (or  charge  conjugation  symmetry),  the  symmetric  state  |u)  has 
exactly  zero  electric  current  expectation  value 

=  0  (2-22) 

at  all  times,  by  the  symmetry  of  this  state.  Likewise  the  mean  charge  density  (u|  p  \v) vanishes  in  this 


charge  symmetric  state.  Thus  the  state  \v)  defined  by  (2.13)-(2.18)  in  a  constant,  uniform  electric  field 
background  is  an  exact  self-consistent  solution  of  the  semiclassical  Maxwell  eqs. 

V  •  E  =  (uj  p  lu)^  =  0  (2.23a) 

dE 


V  X  B  -  ^  =  (ulJlu)^  =  0 


(2.23b) 
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with  both  sides  vanishing  identically.  This  is  analogous  to  the  maximally  0(4, 1)  symmetric  and  time 
reversal  invariant  CTBD  state  |n)  which  satisfies  the  semiclassical  Einstein  eqs.  in  de  Sitter  space 


with  a  simple  redefinition  of  A,  since  {v\T°‘^\v)r  =  5\  ,  c.f.  Sec.  IV 


B.  Instability  of  the  Maximally  Symmetric  State:  Electric  Current 

The  existence  of  a  state  of  maximal  symmetry  does  not  imply  that  it  is  the  stable  ground  state  of 
either  the  de  Sitter  or  electric  field  backgrounds.  In  the  electric  field  case  the  imaginary  part  of  the 
effective  action  and  spontaneous  decay  rate  of  the  electric  field  into  particle/anti-particle  pairs  was  hrst 
calculated  by  Schwinger  m-  By  time  reversal  invariance  the  imaginary  part  of  the  effective  action  (which 
changes  sign  under  time  reversal)  corresponding  to  the  symmetric  |u)  state  vanishes,  in  disagreement  with 
Schwinger’s  result.  As  a  precise  coherent  superposition  of  particle  and  anti-particle  pairs  for  all  modes,  the 


time  symmetric  state  defined  by  (2.13)-(2.18)  is  a  very  curious  state  indeed,  corresponding  to  the  rather 


unphysical  boundary  condition  of  each  pair  creation  event  being  exactly  balanced  by  its  time  reversed 
pair  annihilation  event,  these  pairs  having  been  arranged  with  precisely  the  right  phase  relations  to  come 
from  great  distances  at  early  times  in  order  to  effect  just  such  a  cancellation  everywhere  at  all  times. 
While  mathematically  allowed  in  a  time  reversal  invariant  background,  it  would  be  difficult  to  arrange 
such  an  artificial  construction  and  fine  tuning  of  initial  and/or  boundary  conditions  on  the  quantum  state 
of  the  charged  field  with  any  macroscopic  physical  apparatus,  and  certainly  it  would  not  be  produced 
with  a  more  realistic  adiabatic  switching  on  and  off  of  the  electric  field  background  in  either  finite  time 
or  over  a  finite  region  of  space  [25].  Nor  does  the  state  \v)  minimize  the  Hamiltonian  of  the  system  which 


is  time  dependent  in  the  gauge  (2.1),  or  unbounded  from  below  in  the  static  gauge  Aq  =  Ez. 


The  above  physical  considerations  and  Schwinger’s  earlier  result  suggest  that  there  should  be  an 
instability  of  the  time  symmetric  state  to  nearby  states  in  which  the  hne  tuned  cancellation  between 
particle/anti-particle  creation  and  annihilation  events  is  slightly  perturbed.  In  order  to  probe  these 


nearby  states  we  return  to  (2.3),  and  express  its  general  solution  in  the  form 

/k(i)  =  Ak  vx{u)  +  Hk  Vx(u) 

with  the  (strictly  time  independent)  Bogoliubov  coefficients  required  to  obey 

|AkP  -  |Hkp  =  1  for  all  k, 

in  order  for  the  Wronskian  condition 


(2.24) 


(2.25) 


(2.26) 
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to  be  satisfied.  The  Bogoliubov  coefficients  (j4k,i?k)  may  be  regarded  as  specified  by  initial  data  /k(to) 
and  /k(to)  at  t  =  to  according  to 

Ak(t„)  =  i 

Bk(to)  =  i  (aI^a  - 


t=to 

(2.27a) 

t=to 

(2.27b) 

The  quantized  charged  scalar  field  operator  (2.17)  may  just  as  well  be  expressed  in  terms  of  these  general 


mode  functions  (2.24)  as 


x)  =  ^  ^  /k(t)  /*k(t) ' 


— ikx 


(2.28) 


where  upon  setting  the  Fourier  components  of  (2.17)  and  (2.28)  equal,  the  corresponding  Fock  space 
operators  ^k’  I’elated  to  the  previous  ones  by 



or  its  inverse 


“k 


Tlk 

i?k  a4* 


4*  _  R* 

-^k 

—Bk  Ak 


ufi 

^-k 


(2.29) 



(2.30) 


Hence  if  we  define  the  state  |/)  by  the  condition 


«kl/)  =  &(l/)  =  o. 

this  state  contains  a  non-zero  expectation  value 

(/|ak^«kl/)  =  \B^\"  =  (/|b-l^-kl/) 


(2.31) 


(2.32) 


of  V  quanta.  Conversely  the  |z;)  state  contains  a  non-zero  expectation  value  of  /  quanta.  Since  both  the 
|r;)  and  general  |/)  states  are  pure  states,  and  each  can  be  expressed  as  a  coherent,  squeezed  state  with 
respect  to  the  other,  it  is  best  not  to  use  the  term  ‘particles’  for  either  of  these  expectation  values,  nor  can 
one  decide  a  priori  which  among  them  is  the  ‘correct’  vacuum.  This  illustrates  the  fact  that  the  question 
of  which  vacuum  state  to  choose  is  not  limited  to  de  Sitter  space  or  gravitational  backgrounds  only,  but 
is  characteristic  of  QFT  in  time  dependent  and  persistent  classical  background  fields  more  generally. 

The  most  general  state  which  is  both  spatially  homogeneous  and  charge  symmetric  is  the  mixed  state 
with  a  density  matrix  and  a  finite  expectation  value  of  /  quanta  [26j,  which  we  denote  by 

Tr  (al^alppN^  =  A^k  =  Tr  (bilbl^  p/,Ar)  (2.33) 
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FIG.  3.  The  u\v\\'^  integrand  in  the  current  in  (2.21)  or  (2.34)  as  a  function  of  u  for  two  different  values  of  A,  c.f 


(2.5).  The  curve  with  large  oscillations  is  for  A  =  1  while  the  oscillations  are  much  smaller  in  the  curve  for  A  =  5. 


Computing  the  renormalized  mean  value  of  the  electric  current  in  this  charge  symmetric  state  we  find 


Tr(j2p/,Ar)  =  2e 


(27r)3 


{kz  +  eEt) 


|/k(t)|"(l  +  2iVk)- 


2u;k(t) 


=  4e 


(27r)^ 


{kz  +  eEt)[[N^  +  |i?k|"(l  +  2iVk)]  kA(n)P  +  (1  +  2A^k)  Re  [A^B*^,vl{u)]  }  (2.34) 


where  we  have  used  (2.22)  and  ( 2.24 )-( 2.25)  in  arriving  at  the  second  expression.  Charge  asymmetric 


states  or  spatially  inhomogeneous  states  with  lower  symmetry  could  be  considered  as  well.  In  a  general 
state  with  i?k  7^  0  or  A^k  7^  0,  the  charge  conjugation  and  time  reversal  symmetry  of  the  background  is 
broken  and  the  current  Tr  P/^n)  7^  0-  Because  such  states  correspond  to  charged  particle/anti-particle 
excitations  that  are  rapidly  accelerated  to  ultrarelativistic  energies  by  the  background  electric  field,  they 
lead  to  persistent  currents  that  do  not  decay  and  which  destabilize  the  constant  electric  field  background 


through  the  semiclassical  Maxwell  eq.  (2.23b) 


To  see  this  requires  only  a  qualitative  understanding  of  the  integrand  in  (2.34).  The  three  terms 
tt|uAp,  uRe  (u^),  and  ulm(u;^)  appearing  in  the  integrand  of  (2.34)  are  shown  as  functions  of  u  for 


several  values  of  A  in  Figs,  [sjj^  In  Fig.  [^the  saturation  of  the  function  u  |uaP  at  large  u  is  the  result  of 
acceleration  of  charged  scalar  particles  to  ultrarelativistic  energies  by  the  electric  field,  where  they  make 
a  constant  contribution  to  the  current  integrand.  Hence  if  A^k  and/or  |RkP  in  (2.34)  is  non-zero  for  any 


range  of  k,  such  modes  will  make  a  contribution  to  the  current  proportional  to  the  phase  space  volume 
f  d^k  which  can  therefore  give  an  arbitrarily  large  (jz)  at  late  times,  u  ^  1. 
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FIG.  4.  The  MRe(z;;^)  integrand  in  the  current  in  (2.34)  as  a  function  of  u  for  two  different  values  of  A.  The  left 
panel  is  for  A  =  1  while  the  right  panel  is  for  A  =  0.1  chosen  to  accentuate  the  asymmetry  in  u  O  —u. 
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FIG.  5.  The  ulm  (v^)  integrand  in  the  current  in  (2.34)  as  a  function  of  u  for  the  same  values  of  A  =  1  (left)  and 


A  =  0.1  (right)  as  in  Fig.  Shown  also  is  the  horizontal  line  at  —  exp(— ttA)  around  which  the  average  amplitude 
of  the  oscillations  are  displaced. 


The  oscillatory  terms  in  the  real  and  imaginary  parts  of  uv^  are  shown  in  Figs.  4j|5  The  envelope  of 
the  oscillations  shows  a  saturation  behavior  at  large  |m|  similar  to  Fig.  For  smaller  A  the  oscillations 
are  significantly  offset  from  the  horizontal  axis,  by  ±exp(— ttA),  showing  that  there  will  also  be  a  net 
contribution  to  the  current  from  modes  with  ^4^  ^  0.  Hence  these  contributions  to  (j^)  can  also 

become  arbitrarily  large  if  the  range  of  k  for  which  is  non-zero  is  large. 


An  interesting  special  case  in  which  to  evaluate  (|2.34|)  is  the  adiabatic  vacuum  state  of  initial  data 

1 


/k(io)  = 


y^2uJk{to) 


/k(to)  =  -fwk  - 


^k 

2i^k 


/k 


t=to 


(2.35) 


and  A^k  =  0.  We  denote  this  pure  state  which  matches  the  lowest  order  adiabatic  vacuum  state  at 
the  particular  time  t  =  to  by  \to)-  With  these  initial  conditions  it  is  shown  in  [9]  that  the  Bogoliubov 
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coefficients  are  given  approximately  by 


0{-kz  -  eEto)  +  e{k,  +  eEto) 
B^to)  ~  Sf  e{-k,  -  eEto)  +  Sr*  0{k,  +  eEto) 


where 


Ain  _  A  out  *  _ 

— 


_  Qout  *  _ 

A  A 


A\3 



e  4 


;r  /  2\  2 


r{|-t)  \x)  r{i-f) 

\  —  y  \  A 

A\4  1  /2A4  64 


r(l-f)  Va;  r(A-f) 


(2.36a) 
(2.36b) 

exp(f-f-f)  (2.37a) 
exp(^-f-^).  (2.37b) 


In  fact,  the  step  functions  in  (2.36)  are  smooth  functions  which  interpolate  between  the  two  limits,  but 


this  simple  approximation  is  sufficient  to  illustrate  the  main  features  of  the  current  expectation  value 
which  is  its  linear  growth  in  time  from  the  initial  time  to- 


Substituting  (2.36)  into  (2.34)  with  N\^  =  0,  changing  variables  from  {kz,k±)  to  (u.  A),  and  making 
use  of  the  fact  that  the  functions  u  |uA(n)p  and  u  Re  [u|(m)]  are  odd  functions  of  u  {c.f.  Figs.  [^-|^,  while 
ulm[u^(u)]  is  even  {c.f.  Fig.  [^,  we  obtain 
46  r°°  ( 


{to\jz\to) 


«  ~  (27r)2 


k±dk± 


'0 


dk,{k,  +  eEt)  \BT\^\vx{u)\^  +  Re^-Rr 


in  T^in*^.2  / 


+  /  dk,{k,  +  eEt)  |Rn>A(u)|2  ^  ji^lAT^Br^*vl{u)] 

J —eEto  ^ 

2  p  poo  p\^2eE  (t— to) 

'  dXlm[AfBt^*]  /  duulm[vliu)] 

Jo 


f 


^  Jm‘^l2eE 

1 


27r2 

lm?l2eE 


f 


\  /O'! 

dXe  ^  /  du[uf[\u)f[\u)] 

Jo 


since  from  (2.37),  or  eq.  (4.20b)  of  ref.  [9]  and  (2.13) 


\Qin\2  |^o'Rt|2 

I  A  I  I  A  I 


out  Tjout*] 


ReiA^ffBflf*]  =  Re[Aff^Bf 


Im[4"Rf  *]  =  -lm[Aff^Bff^*]  =  -AlmRf*  =  ^e 


1  ttA 


Im[uf(u)]  =  — 



(2.38) 

(2.39a) 

(2.39b) 

(2.39c) 

(2.39d) 


A  w  A 

Because  of  the  offset  from  the  n-axis  of  uf^\u)f^\u)  =  — ulm[u^(u)]  by  around  which  the 

oscillations  average  to  zero  {c.f.  Fig.  [^,  for  large  t  —  to  ^  oo  the  u  integral  in  (2.38)  is 


'•V2eE  (t—to) 


du[ufl°\u)f[^\u)] 


/2eE  (t—to) 


due  =  x/2eE  {t  —  to)  e 


(2.40) 
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and  hence  (2.38)  gives  for  late  times 


{to\jz\to)i 


47r3 



(2.41) 


which  is  the  same  result  as  that  of  Eq.  (5.22)  in  Ref.  [9],  which  was  obtained  much  more  naturally  in  the 
adiabatic  particle  basis  by  consideration  of  particle  creation  events.  That  treatment  makes  it  clear  that 
the  growth  of  the  current  is  a  cumulative  effect  of  particle  creation  from  the  quantum  ‘vacuum’  which 
continues  unabated  as  long  as  the  constant  electric  field  is  maintained. 

Thus  there  are  states  for  which  the  current  grows  linearly  with  time  related  to  the  steady  rate  of 
particle  creation  in  a  constant  electric  field  background.  Moreover  it  is  clear  from  the  penultimate  line  of 


(2.38)  that  any  perturbation  of  the  symmetric  |u)  state  with  Bogoliubov  coefficients  of  the  form 

(2.36)  obeying  the  conditions  =  A^^*,  =  B^^^*  of  (2.37),  having  constant  but  non-zero  support 


for  arbitrarily  large  and  negative  will  produce  a  cumulative  effect  on  the  current  similar  to  (2.41),  so 
that  {jz)  continues  to  grow  linearly  with  time  for  arbitrarily  long  times.  This  linear  growth  with  time 
implies  that  however  small  the  coupling  e  and  the  coefficient  \m\A^J^B'^^*]  (which  can  be  enhanced  by 
taking  >  0);  the  current  must  eventually  influence  the  background  field  through  the  semiclassical 


Maxwell  eq.  (2.23b).  Thus  the  symmetric  |u)  state  in  a  fixed  constant  uniform  electric  field  background 


is  unstable  to  perturbations  of  the  kind  (2.36).  If  B\^  has  non-zero  support  up  to  some  large  but  finite 


negative  value  {kz)min  =  —Kz  the  linear  growth  in  (2.41)  will  be  cut  off  at  t  —  to  =  EzjeE  but  still  be 


large  and  produce  a  large  backreaction  through  (2.23b). 

If  one  goes  beyond  the  simple  mean  field  approximation  considered  here,  it  is  also  clear  on  physical 
grounds  that  the  introduction  of  a  single  electrically  charged  particle  into  the  |u)  state  will  cause  it  to 
be  accelerated  by  the  electric  field  to  arbitrarily  large  energies,  which  would  allow  it  to  emit  photons  and 
produce  additional  charged  pairs  resulting  in  an  electromagnetic  avalanche.  Allowing  these  additional 
channels  opened  up  by  self- interactions  makes  the  physical  instability  of  the  symmetric  |u)  state  to  small 
perturbations  more  obvious,  although  that  instability  already  exists  even  without  self-interactions,  in  the 


mean  field  approximation,  as  (2.41)  and  (2.23b)  show. 

This  example  of  the  quantum  states  in  a  constant,  uniform  external  electric  field  shows  quite  clearly 
that  the  most  symmetric  state,  with  the  full  symmetry  group  of  the  background  need  not  be  the  stable 
ground  state  of  the  system.  In  this  case  it  is  well  known  that  the  background  is  unstable  to  particle 
creation.  In  the  accompanying  paper  [9j  we  have  shown  how  the  same  conclusion  follows  in  de  Sitter 
space,  for  essentially  the  same  reasons.  The  treatment  above  shows  that  one  need  not  be  committed  to 
any  definition  of  particles  to  discover  the  instability  of  the  electric  field  background  by  perturbations  of 
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the  symmetric  |r;)  state  which  have  support  at  large  canonical  momentum  \kz\-  For  <  0  this  may 
correspond  to  small  physical  kinetic  momentum  kz  +  eEto  at  some  early  initial  time  to-  The  unlimited 
growth  of  the  physical  momentum  kz  +  eEt  with  time  for  fixed  canonical  momentum  kz  in  terms  of  which 
the  initial  state  is  specified  is  the  essential  feature,  and  this  feature  is  found  in  gravitational  backgrounds 
such  as  de  Sitter  space  as  well. 


C.  Relation  to  Quantum  Chiral  Anomaly  in  Two  Dimensions 


The  linear  secular  growth  of  the  current  in  a  background  constant  electric  field  can  also  be  understood 
through  the  Schwinger  anomaly  in  1  +  1  dimensions  [29].  For  that  comparison  we  drop  the  d^k_L/(27r)^ 


integral  in  (2.34)  to  reduce  to  1  +  1  dimensions,  and  further  set  the  mass  m  =  0.  We  have  then 

e^E 


(j 


z/2d 


TT 


(t-to) 


(2.42) 


at  late  times.  Since  scalars  are  essentially  the  same  as  fermions  in  1  +  1  dimensions  one  can  use  the 
bosonization  results  IMj  for  fermionic  QED  to  express  the  current  in  the  form 


{f)= 

\/7T 


(2.43) 

where  the  antisymmetric  symbol  in  two  dimensions  and  x  is  a  pseudoscalar  field  whose  derivative  is 
the  chiral  current 



1 


TT 


This  current  has  the  well-known  chiral  anomaly  m 


=  ^nx  =  ^  =  — 

/TT  ZTT  TT 


(2.44) 


(2.45) 


in  a  background  electric  field.  The  second  order  eq.  (2.45)  for  x  with  the  anomaly  source  in  a  constant. 


uniform  field  has  solutions  independent  of  z  of  the  form 

1  eE  2 


Substituting  this  value  of  x  into  the  electric  current  (2.43)  gives 


/  ■  \  6  ,  . 

{jz)2d  =  X  = - (A  -  to) 

-v/vr  TT 


(2.46) 


(2.47) 


which  recovers  (2.42).  Thus  the  linear  secular  growth  of  the  current  with  time  in  the  massless  limit 


is  related  to  the  two-dimensional  chiral  anomaly  and  the  particular  z  independent  solution  (2.46)  to 
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the  pseudoscalar  field  eq.  (2.45).  This  particular  solution  to  (2.45)  is  associated  with  the  spatially 


homogeneous  initial  state  condition  (2.35)  and  state  specified  by  the  mode  functions  (2.24)  and  (2.36) 


It  is  interesting  to  note  that  although  the  anomaly  eq.  (2.45)  is  Lorentz  invariant,  because  it  is  an 


inhomogeneous  eq.,  none  of  its  solutions  are  Lorentz  invariant.  Thus  the  maximal  Poincare  symmetry 


of  the  fixed  electric  field  background  is  necessarily  broken  by  the  solutions  to  the  anomaly  eq.  (2.45), 


which  leads  to  a  spontaneous  breaking  of  symmetry  of  the  background,  at  least  in  the  semiclassical 
approximation  and  neglecting  backreaction.  This  may  be  seen  also  from  the  effective  action  corresponding 
to  the  2D  chiral  anomaly  |29l  [32] ,  viz. 


Slnombc]  =  ^  J  J  (fx  =  l  J  (fx 


-XOx+ 


TT 


/ii/ 


(2.48) 


where  \F~^{x,x')  is  the  Green’s  function  inverse  of  the  scalar  wave  operator  □  in  two  dimensions.  As 
is  well-known,  the  usual  construction  of  the  Feynman  Green’s  function  for  a  massless  scalar  in  two 
dimensions  is  infrared  divergent  due  to  the  constant  /c  =  0  mode,  and  consequently  no  Lorentz  invariant 
Feynman  function  exists  in  this  case.  Green’s  functions  x')  obeying  different  boundary  conditions 

exist,  but  these  necessarily  break  some  of  the  continuous  or  discrete  symmetries  of  the  background.  Thus 


the  form  of  the  effective  action  of  the  2D  chiral  anomaly  (2.48),  together  with  the  absence  of  a  Lorentz 


invariant  Feynman  Green’s  function  \F~^{x,x')  due  to  infrared  divergences  is  sufficient  to  conclude  that 
the  maximally  symmetric  state  in  a  uniform  constant  background  =  E  is  sensitive  to  non¬ 

invariant  initial  and/or  boundary  conditions  which  break  that  maximal  symmetry.  The  linear  growth 


of  the  current  found  in  (2.41)  and  reproduced  by  the  solution  (2.46)  in  (2.47)  is  symptomatic  of  that 


necessary  breaking  of  the  maximal  symmetry  of  the  classical  background  by  the  quantum  chiral  anomaly. 

It  is  also  interesting  that  this  connection  with  the  anomaly  of  massless  fields  in  two  dimensions 
survives  in  four  dimensions  and  even  if  the  field  has  a  non-zero  mass  m,  whose  main  effect  is  to  suppress 
the  coefficient  of  the  linear  growth  by  the  Schwinger  tunneling  factor  exp(— 7rm^/ei7).  We  shall  see 
there  is  also  an  interesting  connection  to  a  quantum  anomaly  of  massless  fields  in  four  dimensional  de 
Sitter  space,  a  local  condensate  bilinear  of  the  underlying  quantum  field(s)  analogous  to  X)  and  simple 


arguments  analogous  to  ( 2.43 )-( 2.47)  which  lead  directly  to  the  analogous  conclusion  of  instability  of  the 


symmetric  state  and  breaking  of  maximal  de  Sitter  invariance  in  that  case  as  well. 


III.  0(4)  INVARIANT  STATES  IN  DE  SITTER  SPACE 


Turning  to  our  primary  topic  of  de  Sitter  space,  we  develop  the  quantization  and  discussion  of  possible 
‘vacuum’  states  in  de  Sitter  space  analogously  to  the  electric  field  case  of  the  previous  section.  For  an 
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uncharged  scalar  field  $  satisfying  the  free  wave  equation 


(-□  +m2)$  = 


1  d  ( 


dx°‘  \ 



dx^ ) 


+  M^ 


$  =  0 . 


in  a  gravitational  background,  with  ^  the  curvature  coupling.  The  effective  mass 

=  m^  +  ^R  =  m^  +  12 


(3.1) 


(3.2) 


is  a  constant  since  the  Ricci  scalar  R  =  12i7^  is  a  constant  in  de  Sitter  spacetime.  In  the  geodesically 


complete  coordinates  (1.2)  the  wave  eq.  (3.1)  may  be  separated  into  a  complete  basis  of  functions  of 
cosmological  time  yk{u)  times  Ykimii^),  the  spherical  harmonics  on  S^.  A  unit  vector  on  is  denoted 
by  N  with  coordinates 


N{x,  0,  (p)  =  (sin  X  n,  cos  x)  =  (cos  Xi  sin  x  cos  9,  sin  x  sin  6  sin  cp,  sin  x  sin  9  cos  (p) . 
The  Ykim{Y)  harmonics  are  eigenfunctions  of  the  scalar  Laplacian  on  the  unit  satisfying 


(3.3) 


^3  Yklrrii 


1 


d  .  2  d  I  d  .  ^  d  1 

—  sin  +  ^-^^sm6»— + 

oX  i^X  sm  9  o9  o9  sm  9  0(p^ 


Yklnrii  —  1)  Yklmi  (3-4) 


sm  X  L  ^ 

with  the  range  of  the  integer  A:  =  1,  2, . . .  taken  to  be  strictly  positive,  conforming  to  the  notation  of  |13) 
and  [33].  These  harmonics  are  given  in  terms  of  Gegenbauer  functions  C'^^^_^(cosx)  and  the  familiar 
spherical  harmonics  Yi^i  (h)  in  the  form  |27j 


Yklm,{N)=2U\ 


i l2,k{k  —  I  —  ly.  _ 


vr  (A:  +  /)! 

with  Z  =  0, 1, ...  A:  —  1  and  rrii  =  —I, normalized  so  that 


(sinx)  C'^t;_i(cosx)k)m,(h) 


r-27r 


d'^^k'Vm'^klmi=  dxsin^X  d9  sin  9  dp)YZ,,  ,Ykimi  =  dyk^vidm', 


JS-i  ‘  JO  JO 

Note  also  that  Y^i^^iN)  =  Yki-mi{N). 

The  time  dependent  functions  yk{u)  satisfy 


i,mi 


(3.5) 


(3.6) 


+  3tanhtt  +  (A:^  —  1)  sech^u  +  (7^  +  |) 
du‘^  du  \ 


yfc  =  0, 


where  the  dimensionless  parameter  7  is  defined  by 


7 


M2  9 


m  4  “  ■ 


(3.7) 


(3.8) 


In  the  massive  case  M'^  >  (the  principal  series)  7  is  real  and  positive.  With  the  change  of  variables 


to  z  =  (1  —  zsinhtt)/2,  the  mode  eq.  (3.7)  can  be  recast  in  the  form  of  the  hypergeometric  equation.  The 


fundamental  complex  solution  yk  — )■  Vk'y{u)  may  be  taken  to  be 


Vk^{u)  =  Hck^  (sechu)^"'"^  (1  “  *  sinhu)*^  T  |  — h  ^7, - Z7;  A:  +  1; 


1  —  i  sinh  u 


(3.9) 
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where  F  =  2-^1  is  the  Gauss  hypergeometric  function  and 



1 


r  (A:  +  5  +  iy)  r  (/c  +  2  -  *7) 


(3.10) 


is  a  real  normalization  constant,  fixed  so  that  Vk'y  satisfies  the  Wronskian  condition 


iHa^{u) 



=  1 


(3.11) 


for  all  k,  where  a{u)  =  H  ^  coshu  is  the  scale  factor  in  coordinates  (1.2).  Note  that  under  time  reversal 


u  —u  the  mode  function  (3.9)  goes  to  its  complex  conjugate 


Vk^{-u)  =  vl^{u) 


(3.12) 


for  all  >  0. 


If  0  <  (3.8)-(3.12)  continue  to  hold  by  analytic  continuation  to  pure  imaginary  7  =  ii/, 


with  Vk'y  — >■  Vk,u-  The  mode  functions  (3.9)  reduce  to  elementary  functions  in  the  massless,  conformally 
coupled  case 



2  ■  - 


H 


sechu  (sechu  —  i  tanh  u)^  =  _  cos  7  e 


—ikr) 



(3.13) 


and  in  the  massless,  minimally  coupled  case 

m  =  0,^  =  0,z^  =  |: 


Vu  3  = 
/c  ,2 


H(ksec]iu  +  itanhu)  ,  ,  .  ,  xa-  H(k  cos  ri  +  i  sin  ri)  „ 

^ ^  (sechn  -  i  tanh  uF  =  '  ,  k  =  2,3,... 

yj2k{k^  -  1)  V2fc(fc2  -  1) 


(3.14) 


where  the  conformal  time  variable  rj  is  given  by,  c.f.  Fig. 

r/ =  sin“^(tanhtt)  G  (— 2,  f )  ,  cosr/  =  sechn,  tan r/ =  sinh  u  . 


(3.15) 


The  complex  positive  frequency  modes  Vk^y  of  (3.9),  (3.14)  are  undefined  for  the  case  v  =  k  =  1  since 


the  solutions  of  (3.7)  are  non-oscillatory  in  this  case,  and  must  be  treated  separately  [laiM].  This  leads 
to  the  non-existence  of  a  de  Sitter  invariant  vacuum  state  or  Feynman  Green’s  function  \F~^{x,x')  for  a 
massless,  minimally  coupled  scalar  in  de  Sitter  space  IMIES],  that  is  similar  to  that  for  a  massless  scalar 
in  two  dimensional  flat  space  discussed  in  Sec.  |II  C 


The  scalar  field  operator  $  can  be  expressed  as  a  sum  over  the  fundamental  solutions 

00  k—1  I 

^{u,N)  =  Y,Y1  E  + 

k=l  1=0  mi=—l 


(3.16) 
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with  the  Fock  space  operator  coefficients  satisfying  the  commutation  relations 


V  I)  t 

^klmi  1 


^kk'^ll'^mimi  ■ 


With  (3.6),  (3.11),  and  (3.17)  the  canonical  equal  time  field  commutation  relation 



(3.17) 


(3.18) 


is  satisfied,  where  II  =  =  is  the  field  momentum  operator  conjugate  to  the  overdot 

denotes  the  time  derivative  H  d/du  and  6^{N,  N')  denotes  the  delta  function  on  the  unit  with  respect 
to  the  canonical  round  metric 

The  Chernikov-Tagirov  or  Bunch-Davies  (CTBD)  state  |u)  [2l[3l[6]  is  defined  by 


^klmi 


|u)  =  0  V  k,  I,  mi, 


(3.19) 


and  is  invariant  under  the  full  0(4, 1)  isometry  group  of  the  complete  de  Sitter  manifold,  including  under 
the  discrete  inversion  symmetry  of  all  coordinates  in  the  embedding  space,  — ?•  —X^  {c.f.  Fig.  [^,  or 

{u,N)  — )>  {—u,  —N),  which  is  not  continuously  connected  to  the  identity.  The  Feynman  Green’s  function 
in  this  maximally  symmetric  state  is  invariant  under  0(4, 1)  and  also  coincides  with  that  obtained  by 
analytic  continuation  from  the  Euclidean  for  >  0  with  full  0(5)  symmetry  [5].  As  in  the  electric 
field  example  of  Sec.  |^the  existence  or  construction  of  a  maximally  symmetric  0(4,1)  invariant  state 
does  not  imply  that  this  state  is  a  stable  vacuum. 

Alternative  Fock  representations  in  real  u  time  are  clearly  possible.  For  example,  since  the  general 


solution  of  (3.7)  may  be  written  as  the  linear  combination 


ykiu)  =  Ak  Vkjiu)  +  Bk  v*k Ju) , 


and  normalized  by  (3.11|)  in  the  same  way  by  requiring 
iHa^{u) 


d 


d 



=  iH 


f*  ^  f  _  r  ^  f* 

Jk^Jk  Jk^Jk 


=  \Ak\^  -\Bk\‘^  =  1. 


(3.20) 


(3.21) 


The  general  functions  yk  ^  a  ^  fk  may  just  as  well  be  chosen  as  a  basis  of  quantization  of  the  field  by 

OO  fc— 1  I 


ot 


^u,N)  =  Y,Y.  E  [4r..yfc(^)^^^™,(iV)  +  a”^^y^(u)nl^,(iV) 

k=l  1=0  m=—l 

with  the  Bogoliubov  transformation  between  the  corresponding  Fock  space  operators 


(3.22) 


a 


^klmi 
kl  —mi 


Ak  Bl 
Bk  Al 


a 


1^ 

^klmi 

/t 

kl  —mi 


(3.23) 
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or  its  inverse 


^klmi 

,/t 


A*  —  R* 
-Bk  Ak 


klmi 


kl  —mi 


(3.24) 


^kl  —mi 

analogous  to  ( 2.24 )-( 2.30)  of  Sec. 

3 

The  mode  function  fk  =  a^Uk  also  satisfies  the  eq.  of  an  harmonic  oscillator 

d^fk 


dv? 


+  [(^^  -  i)  sech\  +  7^]  /fc  =  0 


(3.25) 


analogous  to  (2.6),  and  (3.25)  which  is  the  starting  point  for  an  adiabatic  or  WKB  analysis  of  particle 


creation  in  [9j.  Here  we  note  that  because  of  (3.21)  the  commutation  relations  (3.17)  are  also  satished 


f  f  + 

i^kim  1  canonical  field  commutation  relation  (3.18).  Hence  we  may  define  a  state  |/) 


corresponding  to  the  general  solution  (3.20)  of  (3.7)  or  (3.25)  by 

4/m,  l/)  =  0  ^  k,l,mi. 


(3.26) 


for  any  set  of  complex  coefficients  satisfying  (3.21).  Since  the  solutions  yk{u)Ykimi{N)  at  hxed 

k  form  an  irreducible  representation  of  the  group  0(4)  for  any  {Ak,  Bk},  these  states  are  invariant  under 
0(4)  rotations  of  but  not  the  full  0(4, 1)  de  Sitter  group  (unless  =  1,  =  0  for  all  k). 

The  0(4)  invariant  states  are  associated  with  a  preferred  u  time  slicing  which  breaks  the  0(4, 1) 
symmetry.  The  Bogoliubov  coefficients  {Ak,Bk}  and  hence  the  particular  |/)  state  may  be  regarded  as 
specified  by  initial  data  yk{uo)  and  yk{uo)  on  the  u  =  uq  Cauchy  surface  according  to 


d 


d 


Ak{uo)  =  iHa  (  ^l^^Vk  -  Vk^vl 


■du^ 


Bk{uo)  =  iHa^  ~ 


U=UQ 


U=UQ 


(3.27a) 

(3.27b) 


States  with  lower  symmetry  than  0(4)  may  be  obtained  by  considering  Bogoliubov  transformations  more 


general  than  (3.24),  mixing  and  of  different  {klmi).  For  example  if  the  relation  (3.24)  is  generalized 
to  akirrii  =  Ofc//  -m  Bogoliubov  coefficients  are  (non-diagonal)  matrices  in 


k,k'  (but  still  diagonal  in  I, mi),  the  corresponding  states  (3.26)  are  0(3)  invariant  only.  These  are 
appropriate  for  the  static  coordinates  of  de  Sitter  space.  All  states  related  to  |n)  by  exact  Bogoliubov 
transformations  of  this  kind  are  pure  states  and  related  to  each  other  by  a  unitary  transformation  [8], 
whether  they  involve  different  (klmi)  or  not. 


The  expectation  value  of  is  non- vanishing  in  the  general  |/)  state  defined  by  (3.20),  (3.24) 


and  (3.26),  viz. 


»&„.!/)  =  \Bkf 


4’, =  AtBl  =  (/|o"L  l/> 


(3.28a) 

(3.28b) 
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Hence  the  general  |/)  ‘vacuum’  state  apparently  contains  ‘particles’  defined  with  respect  to  the  de  Sitter 
invariant  state  |t;).  However,  the  converse  is  also  true  as  the  expectation  values 


=  \Bk\ 


klmi  ^kl  —mi  1^/ 


(3.29a) 

(3.29b) 


are  also  non-zero,  so  that  the  de  Sitter  invariant  ‘vacuum’  may  equally  well  be  said  to  contain  ‘particles’ 
with  respect  to  the  general  0(4)  basis  states  |/).  Since  each  of  the  |/)  states  in  either  case  is  in  fact 
a  coherent,  squeezed  pure  state  with  respect  to  the  others,  with  all  exact  quantum  phase  correlations 


maintained,  it  is  better  not  to  attach  the  label  of  ‘particles’  to  either  set  of  expectation  values  (3.28)  or 


(3.29),  or  the  term  ‘vacuum’  to  any  particular  state  in  de  Sitter  space  at  this  point.  As  in  the  electric 


field  case,  mixed  states  which  are  0(4)  invariant  can  be  defined  through  a  density  matrix  /5/,Ar  with 

Tr(P/,A^  4LAimi)  =  (3.30) 


and  Nk  =  0  reducing  to  the  pure  |/)  state  defined  in  (3.26).  A  non-zero  Nk  above  the  arbitrary  0(4) 
invariant  |/)  ‘vacuum’  is  also  best  not  identified  with  any  physical  particle  number.  In  previous  works 
and  in  a  companion  paper  to  this  one  IHlEa],  we  give  a  definition  of  physical  particle  number  in  de  Sitter 
space  based  on  adiabatic  or  slowly  varying  positive  frequency  basis  ^ . 


IV.  ENERGY-MOMENTUM  TENSOR  OF  0(4)  INVARIANT  STATES 

The  behavior  of  perturbations  of  the  CTBD  0(4, 1)  symmetric  state  in  de  Sitter  space  may  be  studied 
through  the  energy- momentum-stress  tensor,  and  the  potential  backreaction  effects  on  the  background 


geometry  through  the  semiclassical  Einstein  eqs.  (1.7),  analogous  to  perturbations  of  the  symmetric  \v 


state  and  backreaction  effects  of  the  electric  current  through  the  semiclassical  Maxwell  eq.  (2.23b). 
The  conserved  energy-momentum-stress  tensor  of  the  free  scalar  field  is 


Tab  =  (V,<I>)(VfeC^)  -  A  /"(Ve4>)(Vd<i>)  + 



.2^2 



Rah  -  ^R  -  +  9abT\ 



(4.1) 


If  the  Heisenberg  field  operator  in  the  general  0(4)  basis  (3.22)  is  substituted  into  this  expression,  and 


(3.30)  is  used,  the  expectation  value  of  Tab  in  the  general  |/)  state  may  be  expressed  as  a  sum  over  modes. 
Since  these  states  are  spatially  homogeneous  and  isotropic,  and  0(4)  invariant,  we  find  that 


Tr(p/,Arr“)  = 
MPf,NT))=STp^^^ 


hJ  =X, 


(4.2a) 

(4.2b) 
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are  the  only  non- vanishing  components  of  the  renormalized  expectation  value  in  coordinates  (1.2).  Since 
the  renormalization  counterterms  are  state  independent,  they  may  be  subtracted  from  the  mode  sum  for 
the  de  Sitter  invariant  state  with  =  0  once  and  for  all.  The  renormalized  expectation  value 

{v\Tab\v) has  been  computed  in  the  CTBD  state  [H].  Because  of  its  de  Sitter  invariance  this  expectation 


value  satisfies  (4.2)  with  =  —By.  Collecting  then  the  remaining  finite  terms  which  differ  from  this 


when  Ak  ^  7^  0  in  the  general  0(4)  invariant  mixed  state,  one  obtains 

~  X!/  Re  [AkB'l.  e^]  -|-  [Nk  +  |RfcP(l  +  2A^a,)]  ef  | 

k=l 
00 

+  2Nk)  Re  [AkBlp^]  +  [Nk  +  \Bk\\l  +  2Nk)]  pf  } 



P 


f,N 


=  Pv  + 


27T^ 

1 

27r2 


where  we  have  defined 


ef  =  \vk\'^  +  2hRe[vlvf,]  -|-  -b  h^)  \vk\'^  =  3p"  -b  2m^\vk\'^  , 


fc=i 


Efc  =  +  2h  VkVk  +  {ujI  +  h^)vl  =  3p^  +  2m^vl , 


,2„.2 


2  — 


2  2 
Wfe  =  ^  +  m  , 


h  =  —  =  H  tanh  u  , 
a 


(4.3a) 

(4.3b) 

(4.4a) 

(4.4b) 

(4.4c) 


an  overdot  denotes  Hd/du,  a  =  H  ^  coshu,  we  have  suppressed  the  7  subscript  and  also  set  ^  =  g  (but 
kept  m  7^  0)  in  order  to  simplify  the  expressions.  This  is  already  sufficiently  general  for  our  purposes,  as 


the  general  case  I  7^  g  adds  no  essentially  new  features.  By  using  the  mode  eq.  (3.7)  satisfied  by  Vk  one 
may  readily  check  that  the  renormalized  stress  tensor  is  covariantly  conserved, 

_de,  „  _  .  H  d 


H- 


'/.JV 


+  3h{ef^N  +Pf,N)  =  ^f,N)  +  ‘^^Pf,N  =  0 


(4.5) 


du  ’  aP  du 

so  that  it  is  sufficient  to  focus  attention  on  the  energy  density  for  the  general  0(4)  invariant  state. 

Since  the  renormalization  subtractions  have  already  been  performed  in  defining  the  finite  py  =  —By  in 


the  0(4, 1)  invariant  state  |t;),  the  additional  state  dependent  mode  sums  in  (4.3)  must  not  give  rise  to 
any  new  UV  divergences.  This  implies  that  the  Bogoliubov  coefficients  Bk  and  numbers  Nk  must  satisfy 

lim  rA;'^|Rfc|l  =  lim  =  lim  Ik'^Nk]  =  0  (4.6) 

7.  L  -I  k^OO  ^  7.  L  -I 


k^oo 


k^oo 


so  that  all  of  the  sums  over  k  for  the  remaining  state  dependent  terms  in  (4.3)  converge.  States  |/)  whose 


Bogoliubov  coefficients  satisfy  (4.6)  in  addition  to  (3.21)  are  UV  allowedoi  UV  finite  O  (A)  invariant  states 


36].  Finiteness  and  conservation  are  clearly  necessary  conditions  for  the  expectation  value  Tr  to 


be  used  as  a  source  for  the  semiclassical  Einstein  equations  (1.7).  These  properties  of  remain  valid 

for  all  UV  finite  states,  including  those  of  lower  symmetry,  provided  only  that  the  Bogoliubov  coefficients 


fall  off  rapidly  enough  at  large  k,  as  in  (4.6). 
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FIG.  6.  The  top  panel  shows  the  coefficient  of  the  term  in  the  energy  density  fc^e^/(27r^)  of  (4.3a),  for  m  =  H 
and  fc  =  10  in  units  of  H^.  The  bottom  left  panel  shows  the  real  part  and  the  bottom  right  panel  the  imaginary  part 
of  the  coefficient  of  the  AkB^.  term  in  the  energy  density,  namely  Ree^/(27r^)  and  — Im£^/(27r^)  respectively 


of  (4.3a)  and  (4.4a)  again  ior  m  =  H  and  k  =  10,  in  the  same  units. 


As  in  the  current  expectation  value  of  Sec.  |II]  we  seek  a  qualitative  understanding  of  the  terms 
contributing  to  the  energy  density  in  (4.3a)  and  (|4.4[).  There  are  three  kinds  of  terms  for  a  given  k 
in  a  general  0(4)  invariant  UV  finite  state,  namely  those  multiplying  the  factors  Re(Afci3p,  and 

Im(Afci?^)  respectively.  These  are  plotted  in  Figs,  [^through  Fig.  11  In  Fig.  |^the  three  summands 
in  (4.3),  namely  A:^e^/(27r^),  Ree^/(27r^)  and  —  Ime^/(27r^)  are  shown  in  units  of  for  the  case 


m  =  H  and  k  =  10.  The  terms  multiplying  the  complex  coefficient  in  (4.3a)  are  oscillatory, 

while  the  function  multiplying  the  real  coefficient  +  |Rfc6(l  +  2A^fc)  is  non-oscillatory.  The  main 
difference  between  the  coefficients  of  the  real  and  the  imaginary  parts  of  A^B^  is  that  the  former  is 
symmetric  about  tt  =  0  while  the  latter  is  antisymmetric.  The  plots  also  show  that  the  maxima  of  the 
two  oscillatory  functions  occur  for  |n|  of  order  one,  while  the  maximum  of  the  third,  non-oscillatory 
function  is  at  the  symmetric  point  u  =  0  and  much  larger  in  magnitude.  In  all  three  cases  the  functions 
fall  off  for  large  values  of  the  time  |tt|  where  the  scale  factor  a{u)  is  large. 
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FIG.  7.  The  panels  at  left  show,  for  m  =  the  summand  Re£^/(27r^)  in  the  energy  density  of  (4.3a)  and 


(4.4a)  multiplied  by  a^,  with  a  =  H~^  coshu  the  scale  factor.  The  panels  at  right  show  — a^/c^Ime^/(27r^).  From 
top  to  bottom  the  values  of  k  are  fc  =  1,  fc  =  10,  and  k  =  100.  The  values  of  Ukj  from  (4.7)  are  0.88,  3.14  and  5.44 
respectively.  The  plots  show  that  the  envelope  of  the  oscillations  is  proportional  to  a~^  for  |u|  >  Ukj,  but  that  the 
behavior  changes  markedly  for  |u|  <  Uk-f- 


Since  the  field  is  massive  one  might  expect  that  at  large  values  of  the  scale  factor  the  contributions  to 
the  energy  density  would  scale  like  a~^.  To  illustrate  the  power  dependence  on  the  scale  factor  we  plot 
in  Figs.  and  the  coefficients  of  the  real  and  imaginary  parts  of  (1  +  2Nk)AkB'^  multiplied  by  for 
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FIG.  8.  The  panels  at  left  show,  for  m  =  1077,  the  summand  Re£:^/(27r^)  in  the  energy  density  of  (4. Sal  and 


(4.4a)  multiplied  by  a^,  with  a  =  H  ^  coshu  the  scale  factor.  The  panels  at  right  show  —a^k'^  Ime^/(27r^).  From 


top  to  bottom  the  values  of  k  are  k  =  1,  k  =  10,  and  k  =  100,  with  the  corresponding  values  of  Uk-y  from  (4.7) 
0.087,  0.88  and  3.00  respectively,  where  the  behavior  changes  markedly.  It  is  also  observed  that  the  envelopes  for 
|u|  >  Uky  also  scale  like  k"^. 


k  =  1, 10, 100  and  m  =  H  and  m  =  lOH  respectively.  We  observe  that  the  oscillations  have  an  envelope 
which  does  scale  like  for  large  \u\  and  large  a{u).  The  envelope  also  scales  like  independently  of 
m,  so  that  if  we  were  to  sum  modes  up  to  a  large  but  finite  K,  we  would  expect  a  K^/a^  behavior 
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characteristic  of  a  non-relativistic  gas.  However  the  rapid  oscillations,  particularly  for  larger  values  of  m 
and  k,  highlight  the  fact  that  these  are  highly  coherent  quantum  states,  and  the  energy  density  is  not 
that  of  quasi-classical  particles  in  any  sense. 

The  a~^  behavior  of  the  envelope  of  the  oscillations  also  does  not  hold  for  small  |u|.  As  shown  in 


detail  by  a  WKB  analysis  of  the  mode  eq.  (3.7)  in  the  accompanying  paper  [9],  the  mode  functions  and 
adiabatic  vacuum  state  change  character  around  the  times  u  =  ^Uk'y  where 


Uk-i  =  In 


k"^  -\  + +  k‘^ 


(4.7) 


The  modes  are  non-relativistic  for  |m|  >  Uk^,  but  relativistic  for  |u|  <  Uk^y-  For  a  conformal  massless  field 
m  =  0,  with  Vk  =  of  (3.13),  of  (4.4a)  vanishes  identically.  This  accounts  for  the  much  smaller 


values  of  the  energy  densities  in  Figs.  [7  and  in  the  central  regions  where  —Uky  <  u  <  +Uk'y,  where 
there  is  no  simple  behavior  of  the  envelope  of  the  quantum  coherent  oscillations.  The  maximum  of  the 
oscillatory  terms  occurs  for  all  values  of  k  and  m  investigated  at  |u|  ~  1  in  the  central  region.  This 
maximum  saturates  at  a  value  of  order  one  in  77^  units  for  large  k  ^  1,  as  shown  in  Fig. 
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FIG.  9.  The  maxima  of  the  oscillations  of  k^Re  e^/(27r^)  are  plotted  for  several  values  of  k.  The  saturation  of  the 
value  of  maximum  at  large  values  of  k  is  apparent. 


The  non-oscillatory  A;^e^/(27r^)  term  in  the  energy  density  is  shown  in  Figs.  10  and  11  for  m  =  H 
and  m  =  lOH  respectively,  for  fc  =  1, 10  and  100.  In  the  left  panels  the  term  is  multiplied  by  a^{u)  and 
in  the  right  panels  the  term  is  multiplied  by  a^{u).  It  is  clear  that  in  all  cases  the  contribution  from  this 



term  is  proportional  to  a  '’at  large  values  of  the  scale  factor  and  is  proportional  to  o  near  u  =  0. 


In  other  words,  it  blueshifts  in  the  contracting  phase  of  de  Sitter  space  (and  redshifts  in  the  expanding 
phase)  as  a  non-relativistic  fluid  for  large  |u|  but  as  a  relativistic  fluid  for  smaller  |u|.  At  |u|  =  Uk^y  given 


by  (4.7),  the  energy  density  transitions  from  non-relativistic  to  relativistic  behavior  and  for  smaller  |u| 
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FIG.  10.  The  panels  on  the  left  show  for  m  =  H,  the  fc^£^/(27r^)  term  in  the  energy  density  of  (4. Sal  and  (4.4b) 
in  units  of  multiplied  by  a  factor  of  a^.  The  panels  on  the  right  show  this  same  term  multiplied  by  a  factor  of 
a^.  From  top  to  bottom  the  value  of  k  which  corresponds  to  each  set  of  plots  is  k  =  1,  k  =  10,  and  k  =  100,  with 
values  of  Uk-y  of  0.88,  3.14  and  5.44  respectively,  where  the  behavior  changes  from  non-relativistic  to  relativistic. 


the  physical  momentum  k/a  dominates  the  mass  term  in  the  mode  eq.  (3.25).  In  the  non-relativistic 
region  u  >  \uk-y\  the  k  dependence  is  k'^.  However  the  maximum  of  the  k'^s^  term  always  occurs  in  the 
relativistic  region  —Uk.y  <  u  <  +Uk^,  where  the  k  dependence  is  k'^,  so  that  this  maximum  value  grows 
unbounded  for  /c  S>  1,  in  contrast  to  the  oscillatory  terms  which  are  bounded  for  large  k:  Fig. 
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FIG.  11.  The  panels  on  the  left  show  for  m  =  lOH,  the  k'^ term  in  the  energy  density  of  (4.3a)  and  (4.4b) 
in  units  of  multiplied  by  a  factor  of  a^.  The  panels  on  the  right  show  this  same  term  multiplied  by  a  factor  of 
a^.  From  top  to  bottom  the  value  of  k  which  corresponds  to  each  set  of  plots  is  k  =  1,  k  =  10,  and  k  =  100.  The 
values  of  Ukj  where  the  behavior  changes  from  non-relativistic  to  relativistic  are  0.087,  0.88  and  3.00  respectively 
for  this  TO. 


For  the  strictly  massless  conformally  invariant  scalar  field,  there  are  no  oscillatory  terms  since 
=  0  identically  for  Vk  =  and  hence  there  are  no  terms  linear  in  in  the  energy  density  or 

pressure  of  a  general  0(4)  invariant  UV  allowed  state.  The  only  contributions  come  instead  from  the 
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terms  quadratic  in  the  perturbation  from  the  de  Sitter  invariant  CTBD  state  |r;).  Substituting 


into  (4.4b)  with  m  =  0  gives 


27r^ 


4 


2=0  2TT'^a^ 


(4.8) 


showing  that  the  relativistic  behavior  observed  in  Figs.  10  11  in  the  relativistic  region  —Uk^y  <  u  <  Uk-y 
holds  for  all  u  in  the  massless  case.  This  result  can  also  be  obtained  by  conformally  transforming  from  flat 
space  to  de  Sitter  space  the  exact  stress  tensor  for  a  conformal  field  in  a  state  other  than  the  Minkowski 
vacuum  [Ij.  Although  this  is  exactly  the  behavior  one  would  expect  for  a  gas  of  relativistic  particles,  we 
emphasize  that  these  are  still  coherent  quantum  excitations  of  the  pure  |/)  ‘vacuum’  state,  in  which  the 


exact  phase  relations  of  ( 3.28 )-( 3.29)  are  maintained 


In  all  cases  the  perturbations  from  the  CTBD  state  |u)  fall  off  in  the  expanding  half  n  >  0  of  de  Sitter 


space  but  grow  in  the  contracting  half  u  <  0.  The  maximum  value  at  the  symmetric  point  from  (4.8)  is 
given  by 


-  fjsfmax 


27r‘^ 


1  °° 

^  [Nk  +  \Bk\\l  +  2Nk)]  -  ^  [^k  +  IBkIH'^-  + 


k=l 


(4.9) 


where  we  have  approximated  the  sum  by  an  integral  valid  for  large  kmax  =  K,  the  maximum  value  of 
k  for  which  \Bk\‘^  and/or  has  support  consistent  with  the  UV  finiteness  conditions  (4.6).  For  this  to 


produce  a  significant  backreaction  on  the  classical  de  Sitter  geometry  through  the  semiclassical  Einstein 
eqs.  it  is  necessary  for  this  to  be  larger  than  the  background  cosmological  energy  density,  i.e, 


SvrGe,  ^  >A  =  3H^ 


or 


GH^ 

Sir 


[Nk  +  \Bk\^{1  +  2Nk)] 


K 


coshu 


>  1. 


(4.10) 


Clearly  no  matter  how  small  GH^  <C  1,  or  the  state  dependent  perturbation  in  brackets  is,  as  long  as 
their  product  is  non-zero  there  is  always  a  large  enough  (but  still  finite  K)  for  which  this  inequality  is 
satisfied  at  the  maximum  at  tt  =  0.  Since  all  finite  k  modes  are  redshifted  in  physical  momentum  as 


A:/a  — )•  0  for  a{u)  — >■  oo,  perturbations  satisfying  (4.10)  at  u  =  0  have  vanishingly  small  energy  densities 


at  early  terms  u  — )•  —  oo.  Hence  for  any  finite  GH^  >  0  there  is  a  large  class  of  0(4)  invariant  but  de 


Sitter  non- invariant  states  satisfying  (4.10),  which  give  rise  to  energy  densities  that  are  large  enough 


to  produce  significant  de  Sitter  non-invariant  backreaction  effects  at  the  symmetric  point  u  =  0,  all  of 
which  have  exponentially  vanishing  de  Sitter  non-invariant  energy  densities  at  times  in  the  infinite  past 


at  I-  of  ‘eternal’  de  Sitter  space.  The  physical  momentum  HK/a  corresponding  to  the  condition  (4.10) 
for  significant  backreaction  is  of  order  y/HMpi  <C  Mpp  far  less  than  the  Planck  scale  Mpp  so  that  the 
semiclassical  approximation  is  still  reliable. 
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To  see  how  the  general  condition  (4.10)  for  large  backreaction  and  de  Sitter  instability  is  realized  in  a 


specific  physical  state,  which  is  the  one  determined  by  adiabatically  switching  on  of  the  background  anal¬ 
ogous  to  switching  on  of  the  electric  field  in  the  inhnite  past  m,  one  can  choose  Bogoliubov  coefficients 
corresponding  to  the  0(4)  invariant  in  state  of  [H]  prepared  at  the  initial  time  uq.  In  the  contracting 
phase  of  de  Sitter  space  u  <  0  this  corresponds  to  choosing  the  mode  functions  according  to  the  initial 
data  at  u  =  Uq 

Ak  =  6(—Uk'Y  —  Uq)  +  9{uq  +  Uk'j)  (4.11a) 

Bk  =  B%9{-Uk^-UQ)  (4.11b) 


with  Uk'y  defined  in  (4.7).  The  step  functions  are  again  simple  approximations  to  the  actual  smooth  but 


rapid  change  at  —Uk^y-  Since  Bk  =  t)  for  Uk'y  >  —uq  the  mode  sums  in  (4.3)  are  cut  off  for  k  >  K^{uq) 
where 

K^{uo)  =  yj 72  sinh^  uq  +  \  ~  ^  (4.12) 

for  |uo|  1-  Also 


41^^  = 


y/2  sinh(7r7) 


_ ikir  7r7 

e  2  e  2 


TDin  _ 


Y^2  sinh(7r7) 


ifcTT  _  TT^ 

e  2  e  2 


(4.13) 


so  that 


Ain  jDin*  _ 

^/C7  ^k'y  — 


i(-)> 


2  sinh(7r7) 


9{-Ukj  -  Uq) 


oscillates  in  k.  Because  the  term  is  bounded  in  k  for  any  state  as  shown  in  Fig. 


11 


(4.14) 

and  its  coefficient 


(4.14)  oscillates  in  k  for  this  in  state,  its  contributions  tend  to  cancel  when  summed  over  k  in  (4.3),  and 
are  negligible  compared  to  the  non-oscillatory  term  in  the  energy  density  for  large  K.y{uQ),  hence  large 
|uo|-  Retaining  only  the  latter  we  then  have  approximately 


£in  ~  e.,;  + 


1 


2tt'^ 

—TT'y 


47r^  sinh(7r7) 


K^{uo) 

k=l 

pK-,{uo)  u3 

'  dk 


a^(u) 



1 


StT^  g27r7  —  1  Y  cosh  U  J 


K. 



(4.15) 


which  agrees  with  eq.  (8.17)  of  [9]  and  (4.9)  for  the  particular  choice  of  \Bk\‘^  from  (4.13)  and  Nk  =  0. 
This  energy  density,  which  has  an  arbitrarily  small  value  for  u  — )•  —00,  is  blueshifted  as  a  relativistic 
fluid  and  by  u  =  0  has  grown  large  enough  to  comparable  to  the  background  de  Sitter  energy  density 
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and  hence  significantly  affect  the  background  de  Sitter  geometry.  It  satisfies  the  inequality  (4.10)  for 
significant  backreaction  if 


Ky{uo) 


Stt 

(UP 


(e27r7  _ 


(4.16) 


or  from  (4.12) 


|uo|  >  In 


2  _i)\ 


(4.17) 


7  \GH^ J 

which  can  always  be  satisfied  for  early  enough  uq,  and  non-zero  GH^  and  7  in  eternal  de  Sitter  space. 


V.  CONFORMAL  ANOMALY,  STRESS  TENSOR  AND  DE  SITTER  SYMMETRY  BREAKING 


As  in  the  electric  field  example,  the  quantum  vacuum  instability  in  de  Sitter  space  is  illuminated  by 
consideration  of  a  quantum  anomaly,  in  this  case  the  conformal  trace  anomaly  of  the  energy-momentum 
tensor  [llEH].  The  non-local  covariant  effective  action  that  gives  the  conformal  anomaly  in  four  dimensions 

is  Ha  EHISl  [39] 


Sn. 


tb]  =  \  J J A4^(x,x') 


bF  +  b' 


E 


nR\ 

^  J 


(5.1) 


This  non-local  effective  action  is  the  analog  of  (|2.48|)  for  the  chiral  anomaly  in  two  dimensions  and  the 
f  (i^XyJ—g  f  d?x'\/—g'Rx  □  (x,  x')R'x  effective  action  for  the  conformal  trace  anomaly  in  two  dimensions 

m-  In  four  dimensions  there  are  two  invariants  E  =  *Rabcd  =  Rabcd.R^'^^'^  —  ^RabR°'^  +  R^  and 

F  =  Cabcd.C^'^'^^  =  Rabcd.R'^^^'^  “  ‘^RabR'^^  +  5"^^  Contributing  to  the  non-local  anomaly  with  corresponding 
dimensionless  coefficients  b  and  b'  proportional  to  h  in  the  notation  of  |38] .  Being  non-local  in  terms  of 


the  curvature  invariants  E  and  E,  the  one-loop  effective  action  (5.1 )  contains  information  about  non-local 


and  global  quantum  effects,  i.e.  sensitivity  to  initial  and/or  boundary  conditions,  through  the  Green’s 
function  inverse  A'^^{x,x')  of  the  conformally  covariant  differential  operator 

2  „  1 


A4  =  +  2R^^VaVb  -^Rn\  +  PV^R)Va  =  Va  (V“V^  +  2R^^  -  iRg^’^)  Vb  ■ 

00  '  / 


(5.2) 


To  (5.1)  it  is  possible  to  add  any  conformally  invariant  action  (non-local  or  local)  which  does  not  affect 


the  anomaly.  However,  only  the  conformal  breaking  (5.1 )  term  in  the  effective  action  needs  be  retained  in 
a  low  energy  classification  of  operators  in  the  effective  action  [TT|,  and  can  have  relevant  infrared  effects. 


Moreover,  the  effective  action  of  the  anomaly  (5.1 )  is  distinguished  by  being  responsible  for  additional 


massless  scalar  degree(s)  of  freedom  in  low  energy  gravity,  not  present  in  the  classical  theory  |41j.  as  seen 
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also  in  two  dimensions  by  the  shifting  of  the  central  charge  from  N  —  26  to  N  —  25 


this  is  made  explicit  by  rewriting  (5.1)  in  the  local  form 

c  —  ciE)  _a_uc{F) 

iJanom  —  ^  ^anom  '  ^‘^anom  i 


In  four  dimensions 


(5.3) 


by  the  introduction  of  at  least  one  additional  scalar  field,  where  for  example 

=  \  I  d^xV^{-{nip)^  +  2(^R‘^’^-lRg-’’){Va^XVbip)  +  {E-lnR)^}  (5.4) 


is  the  term  related  to  E  in  terms  of  the  additional  scalar  field  (p.  This  scalar  (analogous  to  y  of  Sec.  II C ) 
is  a  new  effective  degree  of  freedom,  not  to  be  confused  with  the  original  scalar  field  <h,  which  describes 
two-particle  correlations  or  bilinears  (relativistic  ‘Cooper  pairs’)  of  the  underlying  scalar,  fermion  or 


vector  QFT  [IHl  HI].  QFT’s  of  different  spin  may  all  be  studied  via  the  effective  action  (5.4),  since  the 
only  dependence  upon  spin  for  free  fields  is  through  the  trace  anomaly  coefficients  b'  and  b,  where 


6'  =  - 


1 


(Ns  +  llNp  62Ny) 


(5.5) 


360(47r)2 

is  the  coefficient  for  the  E  term  in  the  conformal  anomaly  for  non-interacting  scalar  (5),  fermion  (T), 
or  vector  {V)  fields  respectively.  The  b  term  in  the  anomaly  proportional  to  E  gives  rise  to  an  effective 


action  sinlm  similar  to  (5.4)  but  is  less  important  in  de  Sitter  space  where  F  =  CabcdC°‘^‘^'^  =  0  (HHEj. 


Formally  solving  (5.6)  for  (/?  in  a  general  metric  by  inverting  A4,  i.e.  finding  its  Green’s  function 
Aj^(x,  x'),  and  substituting  the  solution  into  (5.4)  returns  the  b'  term  of  the  non-local  form  (5.1).  In  de 


Sitter  space  we  also  have  E  =  24Ff^,  \3R  =  0,  and  the  operator  A4  factorizes,  so  that  the  variation  of 


(5.4)  with  respect  to  (p  yields  the  linear  eq.  of  motion, 


^Mds  =  -n(  -  n+2H^)p  =  I  -  ^  =  12H^ 


(5.6) 


with  a  constant  source.  Because  of  this  constant  source,  analogous  to  (2.45),  and  the  fact  that  the  only 


invariant  scalar  in  de  Sitter  space  is  a  constant,  it  is  clear  that  no  de  Sitter  invariant  constant  solution  to 


(5.6)  for  p  exists.  In  the  local  form  of  the  effective  action  (5.4),  the  freedom  to  add  homogeneous  solutions 
to  (5.6)  is  equivalent  to  that  of  specifying  the  particular  Green’s  function  inverse  Af^{x,x')  dependent 


upon  initial/boundary  conditions  in  the  non-local  form  (5.1).  It  is  also  clear  from  the  factorized  form 


(5.6)  of  A4  in  de  Sitter  space  that  its  inverse 

1 


.-11 


\dS 


2H^ 


[(-□)“^  -  (-□  +2i/2)-i] 


(5.7) 


cannot  be  de  Sitter  invariant,  since  it  is  proportional  to  the  difference  of  the  inverses  of  a  massless, 
minimally  coupled  (^  =  0)  scalar  and  a  massless,  conformally  coupled  (^  =  g)  scalar,  and  no  de  Sitter 
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invariant  form  of  the  former  exists  [34],  Thus  the  breaking  de  Sitter  invariance  and  infrared  sensitivity 
to  initial/boundary  conditions  is  already  apparent  from  either  the  non-local  one- loop  effective  action 


5.ip  and  nonexistence  of  a  de  Sitter  invariant  Feynman  Green’s  function  (5.7),  or  equivalently  from  the 


non-invariance  of  the  solutions  to  (5.6)  and  hence  those  of  the  local  effective  action  (5.4). 


The  form  of  the  breaking  of  de  Sitter  invariance  may  be  studied  through  the  stress  tensor  corresponding 


to  the  local  effective  action  (|5.4|),  whose  variation  with  respect  to  the  metric  gives  the  tensor 
2 


^ab  Vp\ 


<55, 


{E) 


dS 


^9' 


ab 


dS 


=  -2(V(„v9)(Vfe)n(/?)  -b2(VcV7)(V''VaVby?)  -h  2  (□</?) (VaV^vp) 


-^V,V4(Vv^)2]  -AH\Va^){Vb9^)  +  \9ab  -(□7^)"  +  ^n  [{V +  2H\V 
-  ^  VaVfc  □  (^  -h  477^  VaVbV?  -  QabUif  +  SH'^Qab 


(5.8) 


which  is  covariantly  conserved  by  the  use  of  (5.6).  In  (5.8)  we  have  evaluated  Eab  in  de  Sitter  space 


and  used  the  notation  (V(^)^  =  g°'^ iy a^)iy b^)  ■  The  stress  tensor  =  b'Eab  evaluated  on  solutions 
ip  satisfying  the  classical  linear  eq.  (|5.6[)  may  be  used  to  evaluate  the  renormalized  expectation  value 
{Tab) R  of  the  underlying  QFT.  This  is  exact  up  to  state  dependent  (but  curvature  independent)  terms  if 
the  spacetime  is  conformally  flat  as  is  de  Sitter  space,  and  the  QFT  is  classically  conformally  invariant 


We  show  below  in  particular  that  (5.8)  reproduces  the  CTBD  state  value  exactly  for  classically 


conformally  invariant  fields  of  any  spin  with  an  appropriate  choice  of  p. 


Since  the  fourth  order  linear  operator  A4  in  (5.6)  factorizes  into  two  second  order  wave  operators  for  a 


conformally  coupled  and  minimally  coupled  massless  scalar  in  de  Sitter  space,  the  general  homogeneous 


solution  of  (5.6)  in  coordinates  (1.2)  is  easily  found  in  terms  of  lYkimi  and  Vj_.  aYkimi  and  their  complex 


conjugates.  Inspection  of  these  solutions,  (|3.13|)-(|3.15|)  shows  that  the  functions  1  and  3  may  also 
be  written  as  a  linear  combination  of  exp[— f(A:  ±  I)??]-  The  reason  this  rearrangement  of  the  solutions  is 
possible  is  a  consequence  of  the  conformal  properties  of  the  operator  A4.  In  de  Sitter  spacetime  (and  in 
fact  any  conformally  flat  spacetime)  there  is  a  second  factorization  of  A4  into  two  second  order  operators. 


reflecting  the  fact  that  for  a  fixed  k  (5.6)  may  also  be  written 


A4 1  Pk  iv)  Yklrai  {N)  =  cos"^  ri 


Ids 


dr]" 


+  {k-lf 


^2 

dr] 


2  -|-  (A:  -|-  1)^ 


V>kiv)YkimiiN)  (5.9) 


in  conformal  time  t],  where  H'^cos^r]  =  a  Thus  the  homogeneous  solutions  of  (5.6)  are  clearly  lin¬ 


ear  combinations  of  exp[— z(/c  ±  l)r]\Ykimi{N)  and  their  complex  conjugates.  To  these  one  must  add  a 


particular  solution  of  the  inhomogeneous  eq.  (5.6),  which  is  easily  found  in  coordinates  (1.2)  to  be 


Pq  =  21n(coshtt) . 


(5.10) 
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This  particular  solution  is  0(4)  invariant  but  not  0(4, 1)  invariant.  Other  choices  correspond  to  states 


of  lower  symmetry,  but  some  choice  must  be  made  since  the  inhomogeneous  term  in  (5.6)  disallows  the 


de  Sitter  invariant  choice  of  constant  if.  Then  we  may  express  the  general  solution  of  (5.6)  in  the  form 

OO  fc  — 1  I 



k=2  1=0  mi=—l  . 


^klmi 


y^2k{k-l) 


e-iik-i)v  Y  + 


bki 


mi 


■\j2k{k  +  1) 


e-iik+i)i^Y,^.+c.c 


(5.11) 


where  fi{u)  is  the  general  solution  of  (5.6)  for  k  =  1,  constant  on  S'^,  given  by 


fi  (u)  =  2  ln(cosh  u)  +  Co  +  ci  sin  ^  (tanh  u)  +  C2  sech^  u  +  C3  tanh  u  sech  u 
=  2  ln(sec  rj)  +  cq  +  cir]  +  C2  cos^  h  +  C3  sin  rj  cos  rj 


(5.12) 


with  the  Ci  arbitrary  constants  multiplying  the  4  homogeneous  solutions  which  are  functions  only  of  u 


or  conformal  time  rj  defined  in  (3.15).  The  normalizations  of  the  k  >  1  solutions  in  (5.11)  are  chosen 
to  correspond  to  a  previous  canonical  analysis  on  the  conformally  related  Einstein  static  cylinder  M  (8) 
where  the  f  =  2a  field  was  quantized  and  the  (akimi ,  bkimi )  obey  canonical  commutation  relations  (the 
bkimiib\i^^  with  negative  metric)  [H].  Here  we  treat  all  the  expansion  coefficients  {ci,akimiibkimi)  of  the 
general  solution  (5.11)  to  (|5.6|)  for  the  effective  action  in  de  Sitter  space  as  c-numbers. 


For  0(4)  invariant  states  the  stress  tensor  can  only  be  a  function  of  u.  Because  of  the  terms  linear 


in  f  in  (|5.8|)  this  corresponds  to  choosing  all  the  coefficients  akimi  =  =  0  in  (5.11)  for  k>  1.  With 


f  =  fi{u)  substituted  into  (5.8)  we  obtain  the  energy  density 

1 


-  EWfi{u)]  =  if  I  fi  -  -  fi+2hfi  fi  +  ^\^2h  + -h^  -  H^jff  -  2hfi  +  2{H^-3h^)  ipi-Qh  (h  +  i^  )  fi 

(5.13) 


where  a  dot  denotes  the  derivative  H  ^d/du.  Substituting  (5.12)  into  this  expression  gives 


2h' 

e  =  -b'E\[f,{u)]  =  -6b' +  _  (c?  -  ci  -  ci  +  4) 


(5.14) 


The  first  term  gives  the  constant  value  of  the  renormalized  =  —6b'H^  for  the  de  Sitter  invariant  state 
of  a  free  conformal  field  of  any  spin,  with  the  corresponding  pressure  =  6b' H‘^.  The  second  a~^ 


term  shows  that  exactly  the  term  corresponding  to  the  relativistic  limit  obtained  in  Sec.  IV  from  detailed 


analysis  of  the  renormalized  expectation  value  of  the  stress  tensor  of  a  quantum  field  in  the  general  0(4) 


invariant  state  is  simply  reproduced  by  the  anomaly  stress  tensor  (5.8)  with  a  classical  effective  field 
f  =  fi{u).  The  spatial  components 


Eij[fi{u)\  =  6H^gij  +  ^  gij  (cj  -  cl -  cl +  4) 


(5.15) 
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and  eq.  of  state  p  =  e/3  for  the  second  term  are  just  that  required  by  covariant  conservation  (4.5)  for 
this  general  0(4)  invariant  state. 


In  (5.14)  the  arbitrary  coefficients  q  of  the  homogenous  solution  in  (5.12)  appear  and  may  be  related  to 


the  sum  over  the  state  dependent  coefficients  N^iBk  and  K  of  (4.9).  The  de  Sitter  invariant  expectation 
value  for  {Tab) ^  is  recovered  iff 


2  ,  A  2  I  2 

Cl  +  4  —  C2  +  Cg 


(de  Sitter  invariant  (Tab) ) 


(5.16) 


so  that  no  relativistic  radiation  a  ^  term  is  present.  Any  solution  of  (5.6)  of  the  form  (5.12)  with  the 


condition  (5.16)  on  the  coefficients  c*  may  be  taken  as  corresponding  to  the  CTBD  state  and  de  Sitter 


invariant  stress  tensor  with  Sy  =  —pv  It  is  interesting  to  note  in  passing  that  for  the  particular  values 
Cl  =  —2i  and  cq  =  C2  =  cs  =  0,  exp[(/?i(u)]  is  just  the  (complex)  conformal  transformation  that  maps 
flat  space  and  its  Minkowski  vacuum  to  de  Sitter  space  and  the  CTBD  invariant  state  |u).  However,  if 


we  restrict  (^i(u)  of  (5.12)  to  be  real,  and  invariant  under  time  reversal  u  •H'  —u,  corresponding  to  the 


discrete  inversion  symmetry  of  the  CTBD  state,  then  ci  =  C3  =  0,  and  from  (5.16)  C2  =  ±2  so  that 


ip{u)  =  2  ln(cosh  u)  +  Co  ±  2  sech^  u  =  —2  ln(cos  t/)  +  cq  ±  2  cos^  rj 


(5.17) 


is  the  background  solution  to  (5.6)  with  =  —6b'H‘^  =  —py  most  closely  corresponding  to  |u).  Since 


the  stress  tensor  (5.8)  depends  only  upon  derivatives  of  tp,  the  constant  cq  is  irrelevant  and  may  be  set 


to  zero,  so  that  the  choice  of  solution  (5.17)  is  determined  up  to  the  sign  of  the  last  term. 


This  (p{u)  in  (5.17)  is  a  kind  of  mean  value  condensate  of  the  ip  effective  field  in  de  Sitter.  Although 


itself  not  de  Sitter  invariant,  it  gives  a  stress  tensor  corresponding  to  the  de  Sitter  invariant  CTBD  state 
of  the  underlying  QFT.  It  seems  that  one  has  to  consider  more  complicated  expectation  values  such  as 


(TabTcd)  in  order  to  see  directly  the  de  Sitter  breaking  effects  of  the  inhomogeneous  solution  to  (5.6). 


This  is  similar  to  the  de  Sitter  invariant  stress  tensor  {Tab)  obtained  for  a  massless,  minimally  coupled 
held  in  de  Sitter  space  despite  the  non-de  Sitter  invariant  vacuum  state  [34]. 

A  small  variation  of  C2  away  from  ±2  produces  a  de  Sitter  non-invariant  stress  tensor  of  the  form 


(5.14)-(5.15)  which  is  inhnitesimally  small  at  asymptotic  past  inhnity  /_  because  of  its  a  ^  dependence 


upon  the  scale  factor,  but  which  grows  to  hnite  values  at  the  symmetric  time  u  =  0.  The  Ci  satisfying 


(5.16)  are  clearly  a  subset  of  a  wider  class  of  a  three  parameter  family  corresponding  to  0(4)  invariant 


but  non-0(4, 1)  invariant  states.  In  this  parameterization  the  condition  (4.10)  that  the  perturbations  of 
the  CTBD  state  produce  a  large  enough  backreaction  at  u  =  0  to  affect  the  classical  geometry  is 


l67TGH^\b'{ci-cl-4  +  4)\  >  1 


(5.18) 
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Clearly  there  are  a  large  class  of  such  states  all  which  all  have  exponentially  vanishing  de  Sitter  non¬ 
invariant  energy  densities  at  times  u  — )•  — oo  in  the  infinite  past.  Since  a  perturbation  of  the  CTBD 


state  with  infinitesimally  small  energy  density  at  I_  with  coefficients  c*  satisfying  (5.18)  produces  a  large 


backreaction  on  the  geometry  at  tt  =  0,  we  conclude  that  the  de  Sitter  invariant  |u)  state  is  unstable  to 
such  state  perturbations  in  the  initial  data  of  eternal  de  Sitter  space. 

Thus  the  anomaly  effective  action  and  stress  tensor  gives  the  same  result  of  instability  of  the  CTBD 
state  to  perturbations,  obtained  previously  for  massive  scalar  fields,  without  any  need  of  renormalization 


subtractions  or  mode  sums,  although  the  connection  to  the  large  K  cutoff  in  (4.9)  or  (4.10),  or  to  particle 


creation  in  the  in  state  of  (4.15)  or  [9|  is  no  longer  transparent  in  (5.18).  The  anomaly  derivation  of 


the  instability  condition  (5.18)  emphasizes  its  generality,  independent  of  the  particular  case  of  a  non¬ 


interacting  scalar  field,  so  that  (5.18)  holds  for  fields  of  any  spin  simply  by  changing  b'  according  to  (5.5), 


or  more  generally  for  interacting  QFT’s  as  well  with  the  appropriate  b' .  This  result  and  the  composite 
effective  field  (p  is  similar  to  the  generality  of  the  axial  anomaly  derivation  of  the  linear  growth  of  the 


current  in  a  persistent  electric  field  background  in  terms  of  the  bosonized  effective  field  x  (2.47) 


VI.  STATES  OF  LOWER  SYMMETRY:  SPATIALLY  INHOMOGENEOUS  STRESS  TENSOR 


The  expansion  (5.11 )  of  the  anomaly  scalar  p  also  enables  a  general  study  of  states  of  lower  than  0(4) 


symmetry  simply  by  allowing  any  of  the  parameters  akim,  or  bkim,  in  the  general  solution  (5.11)  to  be 


different  from  zero.  Substituting  that  general  solution  for  p  in  (5.8)  gives  a  Tab  which  is  a  function  of 


directions  N  on  as  well  as  u.  In  order  to  study  the  effect  of  these  0(4)  breaking  terms,  we  linearize 


the  anomaly  stress  tensor  around  the  solution  p{u)  of  (5.17)  with  a  de  Sitter  invariant  stress  tensor  by 


if  =  ip{u)  +  (l>{u,  N) 


(6.1) 


for  (j)  a  general  solution  of  (5.6)  with  A4<p  =  Oj  expressed  as  the  sum  of  modes  (5.11).  To  first  order  in 

Eab  =  -2(V(a(^)(Vf,)n(/))  -  2  {V (^aE\<p){Vb)(j))  +  2{VcP){VVaVb<P)  +  2  (y^'V aVb<p)(yc<P) 

+  2{a(p){VaVb(j))  +  2{VaVb‘p)ia(j))  -  ^VaVb  c(p)y d(ti)  i^aT)yb)4>) 

+ffa6{-(n<^)(n<())  +  5n  y)y d<l^)\  +  2H^ y)y d(t>)] 


-  i  '^a'^b  □  <(>  +  4i7^  VaVfe0  -  □  0 

which  is  both  covariantly  conserved  and  traceless.  Using 


(6.2) 


=  P 


(6.3a) 
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VrVup  =  Vi{ip-  hip) 

Vt-Vt-Vt-  ip  =  ip 

Vt-Vt-Vj  ip  =  Vi{ip  —  2h(p  —  hip  +  h?ip) 


(6.3b) 

(6.3c) 

(6.3d) 


operating  on  any  scalar  function  in  the  de  Sitter  metric  in  coordinates  (1.2),  and  the  fact  that  ip{u)  is  a 
function  only  of  u,  we  obtain 


=  2^{U(t))'  +  2{U(p)'  ^  +  2^'^+2'ip^-2{U(p)'<j)-2^{U(t))  -  (□(^)(n(/>) 
^  0)  +6i7^<^(^-  +  ^  (□(/>)"- 4^20  -  □(/) 


(6.4) 


for  the  a  =  b  =  u  component  of  this  linearized  stress  tensor.  Since  the  off-diagonal  metric  components, 
gui,  are  zero,  it  is  slightly  easier  to  compute  in  this  case 



(6.5) 


where 

= —ip  —  (dip)  (j)  —  2(p  —  2h  (j)  +  {h?  —  h)  (f)  +  2  [dip)  —  h4>) 

+  ^  {^^(t>  +  ^<i>-  hip^  -  (p(j)  -  ^  [(□</>)'  -  h  d4>]  +  h<j))  (6.6) 

The  energy  density  can  be  obtained  from  this  component  by  using  the  conservation  equation  =  q 

with  or 


together  with  tracelessness 


ff  +  3/!  (eW  +  pW)  =  6' V'l 

ou  \  / 


/  ^3 


„{1)  -  ip(l) 

p  3  t 


(6.7) 


(6.8) 


so  that 


-(1) 


b' 


^klmi 



,  (t"  - 1) 

Ha* 



for  the  linearized  energy  density  perturbation  in  a  given  {k,l,mi)  mode. 


Substituting  (5.11)  and  (5.17)  into  (6.6),  we  obtain  in  particular  the  contributions 
2  .  ..  2H^ 


3 


secmu  tanhu  (1  =F  2secmu)  x 


I  ^  +  c.a. 


+  ... 


(6.9) 


(6.10) 
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and 


2.  2if3  ,2 

--(□</>)  =  — — sech"^n 

O  O 


’■  Yklm,  +  i  ^ 


—I  _  e 
^/2k 


V2k 


Yklnii  “t“  C.C. 


+  . 


(6.11) 


both  of  which  are  leading  in  k,  where  the  ellipsis  and  all  other  terms  in  (6.6)  are  snbleading  in  k  for 


k  ^  1.  Since  these  terms  are  linear  in  k  for  large  k,  and  because  of  the  additional  factor  of  /c^  from  (6.9), 
these  leading  terms  in  k  in  are  proportional  to  k^.  Taking  into  account  the  time  dependence  next, 


we  observe  that  since  =  sechn  ±  i  tanh 


u 


2  (1) 

as  n  — )•  — oo,  the  leading  sech  u  behavior  of  Vi 


klrri] 


cancels  in  the  sum  of  (6.10)  and  (6.11),  so  that  the  integrand  of  (6.9)  vanishes  at  its  lower  limit,  making 


the  integral  convergent.  The  surviving  subleading  term  then  gives  a  contribution  to  (|6.9|)  of 

du  sechn  tanh  u 


-(1) 


p\  / 

^klmi 


yk^V2 

a 

h'y/2 


3 


{Oiklmi  T  Yklnii  T  c.c. 


k^  sech^n 


iS^klmi  T  )  *  YjfiYfii  T  C.C. 


(6.12) 


as  n  — >■  — oo.  The  integral  in  (6.12)  can  be  computed  exactly  for  n  =  0  with  the  result 


r(l) 

'klmi 


-b'H^^k 


u=0 


{p'klmi  “1“  ^klmi)  i  ^klrrti  “1“  C.C. 


(6.13) 


for  fc  S>  1.  Since  the  contribution  of  this  0(4)  breaking  leading  term  in  k  to  the  total  linearized  energy 
density  is 


CXD  k—1  I 

^‘■>  =  EE  E 

k=l  1=0  mi=—l 


r(b 

'klmi 


(6.14) 


it  falls  off  proportional  to  a  °  from  (6.12)  and  hence  a  factor  of  a  ^  faster  than  the  0(4)  symmetric  terms 


in  (5.14)  as  tt  — ?•  — oo.  From  (6.13)  at  n  =  0  its  maximum  value  grows  with  the  maximum  momentum  K 


for  which  the  coefficients  akimi  or  bkimi  are  non-zero  proportionally  to 

.(1) 


u=0 


'll  ^KLmi 

K  k—1  I 

fk=l 


^  ^  k^  ^  -b’H^  I  dkk^  ^  -b'H^  — 

k=l  1=0  mi=—l  ' 


(6.15) 


and  hence  can  easily  exceed  (5.14)  at  the  symmetric  point  n  =  0  if  X  S>  1.  The  backreaction  of  these 


0(4)  breaking  terms  in  the  stress  tensor  becomes  significant  when 

SttGH^  \b'\  —  \aKlm,\  ~  \b'\  —  \bKlm,\  >  1 


(6.16) 


which  is  ever  easier  to  satisfy  for  a  larger  range  of  state  coefficients  as  K  is  increased.  Thus  these  k 
dependent  0(4)  breaking  terms  begin  smaller  and  for  large  enough  K  grow  larger  to  dominate  the  a~^ 


de  Sitter  breaking  0(4)  symmetric  terms  (5.14)  in  the  stress  tensor  as  a  decreases  from  infinity  at 
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We  conclude  that  the  general  0(4)  invariant  ‘vacuum’  states  |/)  defined  by  (3.20)-p.26[)  are  dynami¬ 
cally  unstable  to  producing  large  deviations  in  the  stress  tensor,  even  more  so  than  the  0(4, 1)  invariant 


|u),  the  larger  the  k  of  the  solution  of  (5.11)  considered.  Hence  the  0(4)  symmetry  subgroup  and  spatial 
homogeneity  is  also  spontaneously  broken  in  ‘eternal’  de  Sitter  space.  This  conclusion  which  follows  from 
the  stress  tensor  of  the  anomaly  could  also  be  obtained  by  calculating  the  expectation  value  {Tab)R  of 
the  underlying  QFT  in  0(4)  non-invariant  states  in  de  Sitter  space. 


VII.  CONCLUSIONS 


The  main  conclusion  of  our  analysis  of  possible  states  in  both  de  Sitter  space  and  the  example  of  a 
constant,  uniform  electric  field  is  that  the  most  symmetric  state  in  such  persistent  background  fields  is  not 
the  stable  vacuum  state.  Unlike  fiat  Minkowski  space  where  the  Poincare  invariant  vacuum  is  determined 
by  a  physical  minimization  of  energy,  no  such  conserved  Hamiltonian  bounded  from  below  exists  in  either 
de  Sitter  space  or  in  a  constant,  uniform  electric  field.  Instead  both  of  these  systems  are  characterized 
by  a  mixing  of  particle  and  anti-particle  modes  with  respect  to  any  proposed  Hamiltonian  generator, 
and  are  therefore  unstable  to  spontaneous  particle  pair  creation  from  the  vacuum  BM-  In  each  case  the 
persistent  or  ‘eternal’  background  classical  field  provides  an  inexhaustible  supply  of  energy  to  create  pairs 
at  a  finite  rate  and  subsequently  accelerate  them  to  ultrarelativistic  particle  energies.  In  this  situation 
one  should  expect  the  symmetric  state  to  be  unstable  to  perturbations  and  capable  of  generating  large 
backreaction  effects,  even  in  a  semiclassical  mean  field  approximation  in  which  particle  self-interactions 
are  neglected.  The  study  of  particle  creation  in  real  time  and  the  resulting  vacuum  decay  rate  given  in  an 
accompanying  publication  [U]  is  perhaps  the  clearest  path  to  the  instability  of  ‘eternal’  de  Sitter  space. 

In  this  work  we  have  provided  two  additional  approaches  to  an  analysis  of  the  instability.  These  are 
both  based  not  on  any  particular  definition  of  particles  but  on  the  study  of  perturbations  of  the  maximally 
symmetric  |r;)  states  and  the  conserved  currents  they  produce.  In  the  electric  field  background  this  state  is 
constructed  in  Sec. |l^  and  is  just  as  well  a  self-consistent  solution  of  the  semiclassical  Maxwell  eqs.  (2.23b) 


as  the  0(4,1)  CTBD  state  is  a  self-consistent  solution  of  the  semiclassical  Einstein  eqs.  (1.7),  with  a 
shifted  cosmological  ‘constant.’  In  each  case  there  are  large  classes  of  perturbations  of  the  symmetric 
state  that  produce  an  electric  current  (jz)  or  stress  tensor  {Tab)  that  are  initially  zero  or  negligibly  small, 
but  which  grow  larger  than  any  prescribed  finite  value.  In  each  case  this  is  due  to  the  blueshifting  of 
field  modes  to  ultrarelativistic  energies.  In  de  Sitter  space  this  occurs  clearly  in  the  contracting  phase 
u  <  0,  and  requires  that  backreaction  of  the  energy-momentum  through  the  semiclassical  Einstein  eqs. 
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0  be  taken  into  account.  Hence  the  assumption  of  a  fixed  de  Sitter  background  is  violated,  ‘eternal’ 
de  Sitter  space  is  unstable  to  perturbations  satisfying  (4.10),  which  produce  large  backreaction  effects, 
and  the  classical  0(4, 1)  symmetry  is  broken  by  quantum  fluctuations. 

The  second  approach  to  instability  of  the  maximally  symmetric  state  in  both  the  electric  field  and  de 
Sitter  backgrounds  is  through  the  relation  to  a  quantum  anomaly.  The  chiral  anomaly  and  bosonization 
method  in  the  2D  Schwinger  model  shows  that  the  invariance  of  the  electric  field  background  is  broken  by 


quantum  effects.  In  the  approximation  of  a  fixed  background  field  the  solutions  of  the  anomaly  eq.  (2.45) 


for  the  effective  boson  field  which  are  spatially  homogeneous  predict  the  linear  growth  with  time  (2.47), 
found  also  by  direct  study  of  the  perturbations  of  the  symmetric  state.  Even  at  the  level  of  the  effective 
action  (2.48),  the  appearance  of  the  Green’s  function  of  the  2D  wave  operator  makes  it  clear  that 
there  will  be  infrared  sensitivity  to  boundary  and/or  initial  conditions  associated  with  the  anomaly. 

In  the  gravitational  case  it  is  the  conformal  trace  anomaly  which  produces  long  lived  infrared  effects 
sensitive  to  either  boundary  or  initial  conditions.  It  is  important  that  the  kinematics  of  the  persistent 
de  Sitter  background  will  always  produce  ultrarelatvistic  energies  for  large  enough  k  so  that  the  stress 
tensor  eventually  behaves  like  that  of  a  massless  conformal  field  which  is  described  by  the  stress  tensor 


of  the  anomaly.  From  the  non-local  form  (5.1)  and  the  infrared  properties  of  the  conformal  operator  A4 
and  its  inverse,  it  is  already  clear  without  detailed  calculation  that  de  Sitter  invariance  is  broken.  As 
in  the  2D  chiral  anomaly  one  can  introduce  a  composite  effective  bosonic  field  ip  whose  eq.  of  motion 


(5.6)  has  a  constant  source  and  therefore  possesses  no  de  Sitter  invariant  solutions.  Since  the  eq.  (5.6) 
is  de  Sitter  invariant  but  none  of  its  solutions  are,  the  anomaly  provides  a  mechanism  for  spontaneous 
breaking  of  de  Sitter  symmetry  [45].  The  behavior  of  the  0(4)  symmetric  solutions  which  break  de  Sitter 
invariance  is  easily  found  and  the  same  conclusion  of  the  instability  of  global  de  Sitter  space  to  these  state 
perturbations  follows.  The  anomaly  approach  is  quite  general  and  shows  that  the  same  large  backreaction 
effect  and  instability  to  initial  state  perturbations  occurs  in  de  Sitter  space  for  fields  of  any  spin. 

In  both  cases  it  is  essential  that  moderate  or  small  physical  momenta  are  blueshifted  to  very  large 
physical  momenta,  arbitrarily  large  if  backreaction  is  turned  off  and  the  background  field  persists  indefi¬ 
nitely.  This  unboundedness  and  relation  to  anomalies  is  a  direct  consequence  of  the  infinite  reservoir  of 
arbitrarily  high  momentum  or  short  distance  modes  in  any  vacuum  state  of  QFT  with  no  UV  cutoff.  The 
physical  momentum  which  first  produces  large  backreaction  effects  is  of  order  y/HMpi.  There  is  thus  an 
interesting  interplay  of  UV  and  IR  physics  in  these  effects,  as  has  been  noticed  by  other  authors  [MlllTj. 

The  instability  of  the  de  Sitter  invariant  state  to  large  backreaction  effects  shows  that  the  0(4, 1) 
symmetry  of  global  de  Sitter  space  is  broken  by  quantum  perturbations  of  the  state.  One  can  construct 
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fully  0(4, 1)  invariant  theories  in  eternal  fixed  de  Sitter  spacetime  mathematically  by  continuation  from 
the  Euclidean  order  by  order  in  perturbation  theory  |48j.  By  this  construction  the  very  state  per¬ 
turbations  responsible  for  the  instability  of  de  Sitter  space  in  real  time  are  disallowed  by  the  Euclidean 
regularity  conditions.  If  one  requires  these  regularity  conditions  explicitly  or  implicitly  by  analytic  con¬ 
tinuation  from  and  fixes  the  geometry  to  be  de  Sitter  exactly,  also  disallowing  the  possibility  of 
dynamical  backreaction  through  the  semiclassical  Einstein  eqs.,  it  is  not  surprising  then  to  hnd  no  sign 
of  the  instability  we  have  discussed  in  this  paper,  which  makes  neither  of  those  assumptions.  This  also 
shows  that  it  is  not  matter  self-interactions  per  se  which  are  critical  for  the  instability,  but  rather  the 
boundary  conditions  on  states  and  Green’s  functions  used.  If  Euclidean  boundary  conditions  are  re¬ 
quired,  interactions  lead  to  no  apparent  instability  [l8].  Conversely,  instability  is  seen  immediately  those 
restricted  boundary/initial  conditions  are  relaxed  (even  in  free  QFT)  and  the  backreaction  effects  of  the 
energy-momentum  tensor  on  the  background  de  Sitter  geometry  are  considered. 

The  anomaly  stress  tensor  shows  further  that  there  are  also  spatially  inhomogeneous  perturbations  of 
the  initial  state  which  break  0(4)  symmetry  and  which  vanish  more  rapidly  in  the  inhnite  past  and  grow 
to  larger  values  at  later  times  than  the  0(4)  symmetric  ones.  This  shows  that  there  is  no  0(4)  invariant 
stable  vacuum  state  in  de  global  de  Sitter  space  either.  In  the  particle  creation  language  this  is  quite 
natural  since  localized  particles  and  fluctuations  about  the  mean  are  always  spatially  inhomogeneous.  It 
will  be  necessary  to  study  these  fluctuations  and  spatial  inhomogeneities  in  order  to  describe  accurately 
the  dissipation  of  either  electric  field  or  cosmological  vacuum  energy  into  particle  modes. 


The  spatially  inhomogeneous  perturbations  are  particularly  relevant  for  inflation  and  cosmological 
models  more  generally,  since  they  will  occur  even  in  the  monotonically  expanding  background  of  de  Sitter 


space  in  the  flat  Poincare  sections  (1.4).  Indeed  because  of  the  homogeneity  of  classical  de  Sitter  space 
all  points  on  the  de  Sitter  manifold  are  locally  equivalent,  so  that  particle  creation  and  the  spatially 
inhomogeneous  fluctuations  in  the  geometry  they  cause  are  not  limited  to  the  global  hyperboloid  coor¬ 


dinates  (1.2).  In  fact,  the  same  anomaly  stress  tensor  evaluated  in  the  static  coordinates  (1.5)  of  de 
Sitter  space  shows  that  it  also  describes  sensitivity  to  boundary  conditions  on  the  cosmological  horizon 
r  =  H~^  which  is  entirely  contained  within  the  flat  Poincare  patch  m-  These  new  cosmological  horizon 
modes  associated  with  the  anomaly  scalar  effective  field  and  stress  tensor  have  potentially  large  effects 
on  the  de  Sitter  horizon.  Their  fluctuations  are  capable  of  generating  and  describing  the  anisotropies  in 
the  CMB  |49].  This  suggests  that  spatial  inhomogeneous  models  of  cosmological  dark  energy  within  a 
Hubble  horizon  volume  which  have  only  rotational  0(3)  symmetry  are  relevant  and  may  be  required  for 
a  resolution  of  the  problem  of  cosmological  dark  energy  in  the  present  universe. 
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